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Fluorescence  of  Optical  Materials  Stimulated  by  ArF 
Laser  Radiation  at  193  nm 

DAVID  J.  WREN* 

Department  of  Chemistry,  Kansas  State  University,  Manhattan,  Kansas  66506 


The  fluorescence  spectra  of  four  optical  materials  (Sl-UV 
quartz,  UV  grade  sapphire,  MgFs,  and  BaF>)  have  been  recorded 
using  ArF  laser  radiation  at  193  nm  (6.4  eV)  as  an  excitation 
source.  The  time-resolved  fluorescence  intensities  have  also 
been  measured  at  several  wavelengths  using  a  transient  digi¬ 
tizer  with  2  ns  resolution.  All  of  the  materials  exhibit  a  fast 
fluorescence  which  follows  the  exciting  laser  pulse  and  has  a 
lifetime  r>  <  7  i  2  ns.  This  fast  fluorescence  is  present  as  a 
broad  continum  from  209  to  1000  nm  with  maxima  in  the  UV 
and  the  III.  In  addition,  there  are  fluorescence  components  with 
longer  lifetimes  which  are  characteristic  of  the  individual  ma¬ 
terials.  The  results  suggest  further  opportunities  for  study  of 
optical  materials  using  ArF  lasers  and  indicate  an  area  of  cau¬ 
tion  to  be  exercised  by  experimentalists  employing  ArF  lasers 
for  excitation  experiments. 

Index  Headings:  ArF  laser;  Laser-induced  fluorescence;  Optics; 
Time-resolved  spectroscopy. 


INTRODUCTION 

The  recent  advent  of  the  rare  gas  halide  lasers  has 
provided  experimentalists  with  a  valuable  new  source  of 
intense  UV  radiation.1  With  peak  radiance  powers  on  the 
order  of  10®  W  and  pulse  durations  of  10  ns,  they  are  ideal 
radiation  sources  for  a  variety  of  work  including  photo- 
lytic  reactions,2  multiphoton  ionization,3, 4  two  photon 
excitation  of  high  electronic  states,5  and  fast  kinetics 
initiated  by  absorption  of  one  or  more  photons.  Despite 
the  body  of  work  employing  these  lasers  which  has  been 
published,  there  has  been  no  mention  in  the  literature  of 
the  effects  of  the  RgX  laser  radiation  on  optical  materials. 
For  many  of  the  above  experiments  and  particularly  the 
last,  interference  from  laser-stimulated  fluorescence  in 
optical  materials  can  be  an  important  factor.  In  this 
paper  we  report  fluorescence  spectra  observed  from  four 
optical  materials  stimulated  by  ArF  laser  radiation  at  193 
nm  (6.4  eV). 

The  materials  studied  here  are:  Sl-UV  grade  quartz, 
UV  grade  sapphire,  MgF2,  and  BaF2.  All  of  these  mate¬ 
rials  are  considered  suitable  optical  materials  at  193  nm 
and  are  not  noted  for  fluorescence  in  the  presence  of  UV 
radiation.  However,  all  four  materials  are  observed  to 
fluoresce  with  different  spectra  and  fluorescence  life¬ 
times.  The  fluorescence  observed  from  quartz  is  partic¬ 
ularly  severe.  The  generation  of  fluorescence  by  intense 
radiation  from  an  ArF  laser  is  not  particularly  surprising, 
but  can  pose  problems  for  the  experimentalist  who  wishes 
to  observe  light  from  a  reaction  or  other  process  initiated 
by  an  ArF  laser. 

The  intention  of  this  work  is  to  report  the  fluorescence 
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spectra  observed  for  the  above  materials,  particularly  as 
a  guide  to  other  experimentalists.  The  discussion  of  the 
results  is  limited  in  scope,  and  we  do  not  attempt  to 
identify  the  mechanism  or  states  responsible  for  the 
fluorescence  observed. 

I.  EXPERIMENTAL 

A  schematic  diagram  of  the  experiment  is  shown  in 
Fig.  1.  The  laser  source  was  a  Lumonics  TE-261  rare  gas 
halide  laser  operating  on  the  ArF  193  nm  lasing  transition 
using  a  mixture  of  20%  Ar,  0.4%  F2,  and  79.6%  He  at  a 
total  pressure  of  1500  Torr.  The  gases  were:  Matheson 
prepurified  Ar  (99.998%),  Matheson  technical  grade  F2 
(>97%)  as  10%  F2  in  He,  and  high  purity  He  (99.995%). 
The  laser  was  operated  in  the  static  fill  mode  with  the 
mixtures  prepared  by  backfilling  the  laser  using  a  Tach- 
isto  model  GMS  150  gas  handling  system.  At  44  kV  and 
5  Hz  the  laser  yielded  an  initial  energy  of  -100  mJ/pulse 
with  an  approximately  linear  decrease  in  energy  with  the 
number  of  pulses  to  —50  mJ/pulse  after  —2500  pulses. 
Pulse  to  pulse  stability  was  ±5%.  The  laser  mixture  was 
normally  renewed  after  —2400  pulses.  The  full  width  of 
half  maximum  of  the  laser  profile  was  measured  to  be 
123  cm-1  and  the  laser  pulse  width  was  14  ±  2  ns. 

The  materials  used  in  this  study  were:  (a)  an  Sl-UV 
quartz  flat,  25  mm  diameter,  3  mm  thick,  supplied  by 
Esco  Products;  (b)  a  UV  grade  sapphire  flat  18  mm 
diameter,  1  mm  thick,  supplied  by  Adolf  Meller;  (c)  a 
MgF2  flat,  25  mm  diameter,  1  mm  thick,  supplied  by 
Harshaw;  and  (d)  a  BaF2  lens,  50  mm  diameter,  3  mm 
thick,  {/ 12,  supplied  by  Harshaw.  The  optical  materials 
were  mounted  at  an  angle  of  -30°,  70  cm  along  the  path 
of  the  laser  beam.  At  this  point  the  laser  beam  cross 
section  is  25  mm  by  6  mm.  Specular  reflection  was 
directed  away  from  the  entrance  slit  of  the  monochro¬ 
mator,  and  a  light  shield  prevented  radiation  from  the 
laser  directly  entering  the  monochromator.  Non-specular 
scattering  of  the  laser  light  remains  sufficiently  intense 
to  appear  at  higher  orders  in  the  fluorescence  spectra. 
Fluorescence  from  the  irradiated  materials  was  collected 
by  a  50  mm  f/1  quartz  lens  and  focused  onto  the  entrance 
slit  of  the  monochromator. 

The  fluorescence  spectra  were  scanned  using  a  0.3  m 
McPherson  monochromator  equipped  with  a  grating 
blazed  at  300  nm,  1200  g/mm  for  wavelengths  200  to  600 
nm,  and  a  grating  blazed  at  1000  nm,  1200  g/mm  for 
wavelengths  600  to  1000  nm.  The  light  was  detected  by 
an  RCA  31034  photomultiplier  tube  cooled  to  -20°C  in 
a  Products  for  Research  rf  shielded  housing.  The  current 
generated  by  the  phototube  was  passed  directly  via  1  m 
of  coaxial  cable  to  the  50  input  of  a  Biomation  transient 
digitizer.  The  Biomation  has  a  1024  channel  time  base 
with  a  minimum  sample  interval  of  2  ns  and  a  resolution 
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Fig.  1.  A  schematic  diagram  of  the  experimental  apparatus  for  mea¬ 
surement  of  ArF  laser  induced  fluorescence. 

of  6  bits  for  each  channel.  Triggering  for  the  Biomation 
was  provided  by  an  RCA  1P28  multiplier  which  viewed 
the  laser  beam  directly  via  a  pinhole.  Driven  at  very  low 
voltage  (~200  V)  the  photomultiplier  tube  signal  was 
sufficient  to  trigger  the  Biomation  on  the  leading  edge  of 
the  laser  pulse  and  the  electronics  were  operated  in  a 
delayed  trigger  mode  with  0  delay  time.  Biomation  trig¬ 
ger  jitter  was  ±  6  ns  at  the  fastest  sample  rate.  Both  the 
Biomation  and  the  ArF  laser  were  enclosed  in  grounded 
Faraday  cages  made  of  copper  mesh,  and  the  Biomation 
power  supply  was  isolated  by  a  line  filter.  These  precau¬ 
tions  were  necessary  to  reduce  rf  interference  generated 
by  the  laser  spark  gap  and  gas  discharge.  Random  and 
systematic  rf  noise  were  reduced  to  an  acceptable  level 
but  not  completely  eliminated. 

The  Biomation  was  controlled  by  a  PDP  11/34  com¬ 
puter  which  stored  the  waveform  generated  by  each  pulse 
and  performed  required  data  handling.  In  measuring  the 
time  dependence  of  the  fluorescence  at  a  particular  wave¬ 
length,  typically  200  waveforms  were  acquired  and  aver¬ 
aged,  and  a  background  subtraction  performed  to  remove 
the  rf  noise.  The  fluorescence  spectra  were  recorded  by 
using  the  Biomation  and  computer  as  a  boxcar  integrator. 
For  each  spectral  point  the  computer  acquired  10  wave¬ 
forms  from  the  Biomation,  averaged  them,  subtracted 
the  rf  background  and  integrated  the  result  over  a  pre¬ 
selected  time  interval  of  the  waveform.  The  monochro¬ 
mator  was  scanned  continuously  at  50  nm/min  and  the 
data  points  collected  were  2.5  nm  apart.  The  monochro¬ 
mator  slit  width  was  1  mm  corresponding  to  a  resolution 
of  2.6  nm.  Lack  of  synchronism  between  computer  and 
monochromator,  because  they  were  not  interfaced,  re¬ 
sulted  in  an  error  of +10  nm  after  scanning  400  nm.  Since 
the  fluorescence  spectra  are  primarily  continua,  this  is 
not  a  serious  problem. 

Since  the  laser  was  operated  in  a  static  fill  mode,  the 
intensity  of  the  laser  decayed  during  the  time  required 
for  a  spectrum  to  be  obtained.  Furthermore,  the  laser 
mixture  was  replaced  when  the  gratings  were  changed  at 
600  nm.  The  dependence  of  the  fluorescence  intensity  on 
the  intensity  of  the  laser  radiation  was  measured  for  all 
of  the  samples.  The  method  was  simply  to  follow  the 
integrated  fluorescence  intensity  at  a  particular  wave¬ 
length  as  the  laser  intensity  decayed.  The  result  in  {ill 
cases  found  the  fluorescence  intensity  to  be  first  order  in 
the  laser  intensity.  All  of  the  spectra  which  appear  here 
have  been  corrected  for  the  decay  in  the  laser  intensity 
using  this  result.  The  spectra  have  not  been  corrected  for 


the  wavelength  dependent  response  of  the  detection  sys¬ 
tem. 

II.  RESULTS 

A.  Sl-UV  Quartz.  Sl-UV  quartz  exhibits  the  most 
intense  fluorescence  of  the  four  materials  examined.  To 
the  eye,  the  quartz  flat  gives  off  a  bright  blue-white  light 
when  struck  by  the  ArF  laser  beam.  The  relative  intensity 
of  the  quartz  fluorescence  at  400  nm,  time  integrated 
from  0  to  2000  ns,  is  compared  to  that  of  the  other 
materials  in  Table  I.  The  quartz  emission  is  ~40  times 
stronger  than  the  others,  which  all  have  comparable 
fluorescence  intensities.  The  relative  intensities  are 
wavelength  dependent,  and,  in  general,  the  quartz  is 
always  by  far,  the  strongest.  The  blue  fluorescence  of  the 
quartz  flat  may  explain  the  observation  that  the  ArF 
laser  beam  from  the  Lumonics  laser  is  visible  to  the  eye 
as  bright  blue,  despite  the  fact  that  the  lasing  occurs  in 
the  UV.  The  light  is  likely  due  to  the  fact  that  the  output 
coupler  of  this  laser  is  made  of  quartz. 

The  fluorescence  spectrum  of  Sl-UV  quartz  is  shown 
in  Fig.  2.  The  intensity  recorded  by  each  point  represents 
the  time-integrated  fluorescence  from  0  to  2000  ns.  The 
spectrum  consists  of  a  broad  continum  from  200  to  1000 
nm  with  maxima  at  300  nm  and  at  950  nm.  Several 

TABLE  I.  Relative  fluorescence  intensities. 

Material  Ikkl"'  h 

Sl-UV  quartz  1.0 

Sapphire  0.023  ±  0.005r 

MgF2  0.022  ±  0.004' 

_ BaF2 _ 0.030  ±  0,006* 

‘  Normalized  to  Sl-UV  quartz. 

h  A  —  400  nm,  fluorescence  integrated  from  0  to  2000  ns. 
c  Error  limits  are  ±  0.002  for  the  measured  data.  Differences  in  the 
geometries  of  the  irradiated  samples  could  not  be  accounted  for 
exactly  and  an  overall  error  in  the  relative  intensities  of  20%  is 
estimated. 
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Fig.  2.  Spectrum  of  the  fluorescence  from  an  irradiated  Sl-UV  quartz 
sample.  Fluorescence  is  time  integrated  from  0  to  2000  ns.  Resolution 
is  2.6  nm.  The  narrow  features  at  386,  772,  and  965  nm  are  scattered 
ArF  laser  light  in  second,  fourth,  and  fifth  order.  The  line  in  third  order 
at  579  nm  is  too  weak  to  observe.  Spectrum  is  corrected  for  variation 
in  laser  intensity,  but  not  for  wavelength  response  of  detector. 


628 


Volume  34,  Number  6,  1 980 


narrow  features  are  due  to  scattered  ArF  laser  radiation 
in  second  to  fifth  order.  The  decay  of  the  quartz  fluores¬ 
cence  intensity  at  300  nm  is  plotted  against  time  in  Fig. 
3.  The  curve  exhibits  a  initial  fast  decay  with  a  decay 
time  of  16  ±  2  ns  and  a  second  much  slower  decay  with 
a  decay  time  of  450  ±  40  ns.  The  decay  times  quoted  for 
this  fluorescence  and  all  decay  times  later  quoted  were 
determined  by  logarithmic  least  square  fits  to  the  data 
over  a  range  where  a  logarithmic  plot  of  the  intensity  vs 
time  is  linear.  No  fit  was  made  with  fewer  than  10  data 
points  and  for  long  decays  several  hundred  points  were 
used.  The  error  limits  quoted  are  to  1  standard  deviation 
in  the  slopes  of  the  logarithmic  lines  fit  to  the  data.  The 
fluorescence  decay  measured  at  all  wavelengths  contains 
both  a  fast  (0  to  40  ns)  and  slow  (40  to  2000  ns)  component 
with  variation  in  relative  intensities  of  the  fast  and  slow 
light  and  with  the  decay  lifetime  of  the  slow  emission 
varying  from  173  to  556  ns.  The  slow  emission  present  at 
all  wavelengths  contributes  the  dominant  fraction  of  the 
time  integrated  fluorescence  intensity  and  hence  the 
spectrum  in  Fig.  2  is  primarily  that  of  the  slow  emission. 
No  attempt  was  made  to  separate  the  fast  and  slow 
components  over  0  to  40  ns  in  order  to  obtain  a  spectrum 
of  just  the  fast  emission  component.  The  dependence  of 
the  fluorescence  intensity  on  the  intensity  of  the  laser 
source  was  measured  for  both  the  combined  fast  and 
slow  components  integrated  over  0  to  2000  ns  and  the 
slow  component  alone  integrated  over  40  to  2000  ns.  The 
results  for  both  cases  found  the  fluorescence  to  be  first 
order  in  the  ArF  laser  intensity  independent  of  wave¬ 
length. 

B.  Sapphire.  The  sapphire  fluorescence  spectrum  is 
shown  in  Fig.  4.  It  is  similar  to  that  from  quartz  in  that 
the  spectrum  is  a  continum  from  200  to  1000  nm  with 
maxima  in  the  blue  at  325  nm  and  in  the  red  at  975  nm 
and  also  contains  ArF  laser  lines.  These  lines  appear 
more  intense  than  in  Fig.  2,  due  to  the  lower  intensity  of 
the  sapphire  fluorescence.  The  photomultiplier  sensitiv¬ 
ities  used  to  record  the  spectra  for  sapphire,  BaF2,  and 
MgF2  were  higher  than  that  used  for  the  quartz  spectrum. 
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Fic.  3.  Fluorescence  intensity  at  300  nm  as  a  function  of  time  from 
irradiated  Sl-UV  quartz.  The  decay  time  of  fast  fluorescence  is  14  ±  2 
ns  and  the  decay  time  of  the  slow  fluorescence  is  450  ±  40  ns. 
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Fig.  4.  Spectrum  of  the  fluorescence  from  an  irradiated  UV-grade 
sapphire  sample.  Fluorescence  is  time  integrated  from  0  to  2000  ns.  The 
narrow  features  at  386,  579,  772,  and  965  nm  are  scattered  ArF  laser 
light  in  second,  third,  fourth,  and  fifth  order.  Spectrum  is  corrected  for 
variation  in  the  laser  intensity  but  not  for  wavelength  response  of 
detector. 

As  a  result  the  ArF  line  in  higher  orders  saturates  the 
detection  system  and  appears  off  scale.  The  intensities  of 
the  ArF  lines  appearing  in  the  spectra  of  Figs.  4,  6,  and 
7  are  not  accurate  but  give  a  qualitative  guide  to  the 
relative  intensities  of  the  fluorescence  spectra.  A  unique 
feature  of  the  sapphire  spectrum  is  a  group  of  narrow 
bands  around  700  nm.  These  bands  are  responsible  for 
the  bright  red  fluorescence  of  this  material,  which  is 
visible  to  the  eye. 

Unlike  the  Sl-UV  quartz,  the  sapphire  fluorescence 
spectrum  is  primarily  due  to  a  fast  emission  process.  The 
light  intensity  at  300  nm  plotted  as  a  function  of  time  is 
shown  in  Fig.  5a.  The  fast  emission  follows  the  laser  pulse 
even  more  closely  than  in  quartz  with  a  decay  time  of  7 
±  2  ns.  Similar  decay  curves  are  observed  at  other 
wavelengths  except  those  spanned  by  the  red  bands. 
These  bands  arise  from  a  much  different  fluorescence 
process  and  have  a  very  long  decay  lifetime.  The  fluores¬ 
cence  intensity  as  a  function  of  time  for  emission  at  694.5 
nm  is  shown  in  Fig.  5b.  This  emission  reaches  a  maximum 
0.6  ±  0.2  ms  after  the  laser  pulse  and  decays  with  a 
lifetime  of  1.1  ±  0.2  ms.  The  intensities  of  both  the 
continuum  fast  fluorescence  and  the  slow  fluorescence  of 
the  red  bands  were  found  to  be  first  order  in  the  intensity 
of  the  ArF  laser  intensity. 

C.  MgF2.  To  the  eye,  the  fluorescence  of  MgF2  is  the 
weakest  of  all  the  materials  studied.  There  is  only  a  faint 
blue-purple  emission  visible  when  irradiated  by  the  ArF 
laser.  The  fluorescence  of  MgF2  like  of  that  of  sapphire 
is  dominated  by  a  fast  decay  component  similar  to  that 
shown  in  Fig.  5a.  This  emission  again  follows  the  laser 
pulse  shape  closely  with  a  decay  time  of  7  ±  2  ns  and  is 
present  at  all  wavelengths.  The  fluorescence  spectrum  of 
this  fast  emission  alone  is  shown  in  Fig.  6a  where  the 
intensity  was  integrated  only  over  the  light  emitted  from 
0  to  40  ns.  The  basic  features  of  the  spectrum  are  similar 
to  those  of  sapphire  except  for  an  additional  maximum 
in  the  UV.  The  fluorescence  of  the  band  at  300  nm  is 
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FlG.  5.  a.  Fluorescence  intensity  at  300  nm  as  a  function  of  time  for 
irradiated  sapphire,  b.  Fluorescence  intensity  at  694.5  nm  as  a  function 
of  time  for  irradiated  sapphire.  Note  the  change  in  scale  of  the  time 
axis  compared  to  Fig.  5a. 

complicated  by  the  addition  of  a  second  component  with 
a  decay  lifetime  of  32  ±  2  ns. 

The  fast  emission  from  MgF2  is  much  more  intense 
than  the  slow  emission  and  dominates  a  time-integrated 
spectrum.  However,  there  is  a  weak  slow  decay  fluores¬ 
cence  with  a  lifetime  ranging  from  400  to  495  ns  depend¬ 
ing  on  the  wavelength.  The  fluorescence  spectrum  of  this 
slow  emission  is  shown  in  Fig.  6b  where  a  trigger  delay 
on  the  Biomation  acquisition  has  removed  the  fast  light 
(0  to  40  ns),  and  the  intensity  is  time-integrated  from  40 
to  2040  ns.  The  dominant  feature  of  this  spectrum  is  a 
maximum  at  307.5  nm  corresponding  to  the  shortest 
wavelength  UV  maximum  of  the  spectrum  shown  in  Fig. 
6a.  Also  present  are  very  strong  lines  from  the  scattered 
ArF  laser  line.  These  are  present  here  despite  the  trigger 
delay  because  the  detection  system  has  not  yet  recovered 
from  the  saturation  by  the  intense  signal.  Very  weak 
features  are  present  near  800  nm.  The  band  at  812.5  nm 
has  a  decay  similar  in  shape  to  the  694.5  nm  band  of 
sapphire  but  reaches  a  maximum  in  intensity  at  0.8  ±  0.2 
jus  and  decays  with  a  lifetime  of  2.1  ±5  (is. 
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D.  BaF2.  The  spectrum  for  the  fluorescence  of  BaF2 
is  shown  in  Fig.  7  for  light  intensity  integrated  from  0  to 
2000  ns.  Visibly  the  BaF2  fluoresces  with  a  faint  red  color 
that  is  much  less  intense  than  the  sapphire.  The  decay 
curves  for  the  fluorescence  are  very  similar  to  those  of 
MgF2.  At  all  wavelengths,  except  for  the  continuum  band 
around  300  n,  a  fast  emission  with  decay  time  7  ±  2  ns 
dominates.  The  300  nm  band  derives  from  a  slow  emis¬ 
sion  with  a  decay  time  of  ~30  ns  which  contributes  ~50% 
of  the  intensity  of  the  double  peaked  band. 

III.  DISCUSSION 

The  experimental  results  demonstrate  that  ArF  laser 
radiation  is  particularly  efficient  at  stimulating  fluores¬ 
cence  in  many  and  possibly  all  optical  materials.  The 
case  of  quartz  is  particularly  outstanding.  For  this  ma- 
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Fig.  6.  a.  Spectrum  of  the  fluorescence  from  an  irradiated  MgFj  sam¬ 
ple.  Fluorescence  is  time  integrated  from  0  to  40  ns.  The  narrow 
features  at  386,  772,  and  965  nm  are  scattered  ArF  laser  light  in  second, 
fourth,  and  fifth  order.  Spectrum  corrected  for  variation  in  laser  inten¬ 
sity  but  not  for  wavelength  response  of  detector,  b.  Spectrum  of 
fluorescence  from  an  irradiated  MgF-  sample.  Fluorescence  time  inte¬ 
grated  from  40  to  2040  ns.  Note  the  much  higher  relative  intensity  of 
the  ArF  lines  compared  to  Fig.  6a.  Spectrum  corrected  for  variation  in 
laser  intensity  but  not  for  wavelength  response  of  detector. 
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Fig.  7.  Spectrum  of  fluorescence  from  an  irradiated  BaF_.  sample. 
Fluorescence  time  integrated  from  0  to  2000  ns.  The  narrow  features  at 
386,  772,  and  965  nm  are  scattered  ArF  light  in  second,  fourth,  and  fifth 
order.  Spectrum  corrected  for  variation  in  laser  intensity  hut  not  for 
wavelength  response  of  detector. 

terial  the  blue-white  fluorescence  is  quite  bright  to  the 
naked  eye.  While  it  is  well  known  that  most  fused  silicas 
fluoresce  strongly  in  the  blue  when  exposed  to  UV  light, 
the  purest  grades  of  fused  quartz,  including  the  sample 
used  here,  have  in  the  past  shown  no  fluorescence  when 
exposed  to  conventional  forms  of  UV  light.6  An  exception 
is  the  observation  of  a  strong  emission  band  at  430  nm 
observed  for  both  crystalline  quartz  and  fused  silica  when 
the  samples  are  cooled  to  liquid  nitrogen  temperature 
and  exposed  to  light  at  295  nm.6  The  strong  fluorescence 
observed  here  where  none  is  expected  underlines  the 
efficacy  of  the  ArF  laser. 

The  fluorescence  with  the  slow  decay  lifetime  ~200  to 
500  ns  observed  for  the  irradiated  quartz  is  likely  due  to 
impurities  in  the  material.  Virtually  all  of  the  luminescent 
centers  in  silicate  glasses  have  been  traced  to  the  pres¬ 
ence  of  impurities,  primarily  rare  earth  and  transition 
metals.7  Sigel6  has  suggested  that  the  absence  of  fluores¬ 
cence  previously  noted  for  high  purity  silicas  is  related  to 
the  fact  that  well  known  luminescent  activators  such  as 
copper,  manganese,  and  silver  are  present  in  less  than  50 
ppb.8  At  the  high  power  densities  of  the  rare  gas  halide 
lasers,  however,  such  low  impurity  levels  may  no  longer 
be  sufficient  to  reduce  the  fluorescence.  There  are  high 
purity  grades  of  silica  other  than  that  studied  here,  which 
have  even  lower  metal  contents  (e.g.,  Suprasil-W,  Amer- 
sil,  Inc.)  and  which  may  exhibit  less  fluorescence. 

The  fluorescence  spectrum  of  the  sapphire  sample  is 
composed  of  light  with  a  very  short  lifetime  except  for 
the  narrow  bands  in  the  red.  The  fast  fluorescence  can 
be  traced  to  F  and  F+  centers.  Crystalline  sapphire  is 
known  to  have  absorption  bands  at  6.1  and  6.0  eV  which 
have  been  assigned  to  these  centers.9  10  These  absorption 
bands  are  conveniently  close  to  the  ArF  laser  energy  6.4 
eV.  An  F  center  emission  band  at  413  nm  (3  eV)  has  been 
observed  with  a  lifetime  of  36  ±  4  ns9  which  we  do  not 
observe.  However,  Evans  and  Stapelbroek"'  have  re¬ 
ported  a  326  nm  (3.8  eV)  emission  which  they  assign  to 
the  F*  center  and  for  which  they  measure  a  lifetime  of 


tk  s  7  ns.  This  agrees  very  well  with  the  present  obser¬ 
vation  of  a  fluorescence  band  at  325  nm  with  radiative 
lifetime  tk  =  7  ±  2  ns. 

No  attempt  will  be  made  to  assign  the  fluorescence 
observed  for  the  two  alkaline  earth  fluorides  MgF2  and 
BaF2.  The  fact  that  most  of  the  fluorescence  observed 
for  these  two  materials  has  a  decay  time  equivalent  to 
the  fast  fluorescence  observed  for  sapphire  however  sug¬ 
gests  that  F  center  or  F*  center  transitions  may  be 
responsible. 

The  long  fluorescence  lifetimes  of  the  red  bands  near 
700  nm  in  sapphire  (th  =  1.1  ±  0.2  ms)  with  the  band  at 
812.5  nm  in  MgF2  (th  =  2.1  ±  5  /is)  suggest  that  they  are 
of  very  different  origin  and  possibly  metal  atom  lattice 
impurities.  It  is  interesting  to  note  that  the  maxima  of 
these  emissions  occur  long  after  the  laser  pulse  has  been 
completed  (0.6  ±  .2  ms  and  0.8  ±  0.2  /is,  respectively). 
This  requires  the  existence  of  a  dark  long-lived  interme¬ 
diate  state  or  some  energy  pooling  mechanism. 

A  feature  common  to  all  of  the  spectra  is  a  maximum 
in  the  fluorescence  in  the  IR  near  950  nm.  The  maximum 
is  actually  an  artifact  since  the  spectra  are  not  corrected 
for  wavelength  response.  The  sensitivity  of  the  RCA 
31034  photo-tube  decreases  very  rapidly  near  900  nm  (at 
900  nm  the  response  is  10“'  times  that  at  800  nm),  and 
the  mere  observation  of  a  signal  above  900  nm  indicates 
that  the  actual  fluorescence  is  very  intense  and  may 
extend  far  into  the  IR. 

A  final  observation  is  that  the  fluorescence  intensity 
for  all  four  materials  was  first  order  in  the  ArF  laser 
intensity.  This  means  that  no  multiphoton  absorption 
events  were  responsible  for  the  fluorescence.  At  the  high 
power  density  of  the  ArF  laser,  occurrence  of  a  two 
photon  event  is  not  unreasonable  and,  moreover,  would 
involve  12.9  eV,  energy  sufficient  to  generate  color  cen¬ 
ters  in  all  materials  except  LiF.  The  absence  of  such 
events  noted  here  may  not  hold  for  experiments  where 
the  laser  beams  are  focused  to  provide  even  higher  power 
densities.  Color  centers  generated  by  multiphoton  events 
may  result  in  long-lived  or  permanent  coloration  of  the 
material. 

The  fluorescence  data  reported  here  are  specific  to  the 
optical  samples  studied  and  may  be  expected  to  vary 
depending  on  the  purity  of  the  material.  However,  in  the 
cases  of  sapphire,  MgF2,  and  BaF2  where  most  of  the 
fluorescence  decays  very  rapidly  and  is  probably  due  to 
unavoidable  lattice  defects,  variations  in  sample  purity 
should  not  change  the  overall  fluorescence  spectra.  A 
second  factor  for  consideration  affecting  the  fluorescence 
is  the  past  history  of  the  optical  material.  A  second 
sapphire  flat  which  had  undergone  prolonged  exposure 
to  radiation  from  a  Xe  resonance  lamp  (147  nm)  in  our 
laboratory  was  also  examined.  This  window,  which  had 
no  visible  discoloration  or  damage,  fluoresced  a  very 
intense  blue-violet  with  a  decay  time  of  seconds  when 
exposed  to  ArF  radiation.  The  experimental  results  pre¬ 
sented  here  are  only  a  guide  to  what  may  be  expected  for 
quartz,  sapphire,  BaF2,  and  MgF>. 

IV.  CONCLUSION 

The  preceding  discussion  has  been  very  brief  and 
makes  no  attempt  to  explain  most  of  the  observed  fluo- 
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rescence  features  with  any  rigor.  This  study  was  not  an 
exhaustive  time-dependent  analysis  of  ArF  laser  stimu¬ 
lated  fluorescence;  however,  it  clearly  demonstrates  that 
the  ArF  laser  generates  considerable  fluorescence  in 
most,  if  not  all,  optical  materials.  The  unwary  experi¬ 
menter  is  forewarned  that  this  fluorescence  may  cause 
serious  difficulties,  particularly  in  the  case  of  quartz, 
when  studying  reactions  by  optical  methods  on  time 
scales  of  nanoseconds. 

The  efficiency  of  the  ArF  laser  in  exciting  fluorescence 
in  optical  materials  suggests  that  it  may  be  a  very  useful 
tool  for  studying  the  optical  properties  of  such  materials, 
as  well  as  such  related  properties  as  photoconductivity. 
The  pulse  duration  of  the  ArF  laser  is  short  enough  to 
allow  study  of  lifetimes  for  fluorescence  down  to  the 
nanosecond  region.  Much  more  work  remains  to  be  done 
to  fully  understand  the  events  occurring  when  an  optical 
material  is  irradiated  at  6.4  eV  with  megawatt  power. 
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ABSTRACT 

3 

The  XeCl(B-X)  and  (C-A)  emission  spectra  obtained  from  read  n  of  Xe(  or 
3 

P^)  with  Cl?,  CCl^,  and  COC]^  in  the  presence  of  He,  Ne,  Ar,  Kr  -d  bath 
gases  were  used  to  study  the  vibrational  relaxation  and  transfer  e  ien  the  B 
and  C  states  of  XeCl.  By  using  the  different  Cl  donors  different  ranges  of 
vibrational  energy  were  emphasized.  The  bound-free  emission  spectra  were  simulated 
for  various  pressures  of  bath  gas  to  obtain  vibrational  distributions.  Numerical 
modeling  of  the  XeCl (B)  and  XeCl (C)  vibrational  populations  and  the  B/C  intensity 
ratio  as  a  function  of  pressure  gave  rate  constants  for  vibrational  relaxation 
and  transfer,  as  well  as  the  model  for  the  state-to-state  processes.  For  Ar  as 
the  bath  gas,  vibrational  relaxation  can  be  characterized  by  an  exponential  gap 
model,  P..  a  e  ^  ^T  with  rate  constants  of  1-6,  6-12,  and  20-30  x  10 

ij 

3  -1  -1 

cm  molec  sec  for  the  v  ranges  of  0-30,  30-70,  and  70-130,  respectively.  The 

rate  constants  for  electronic  state  transfer  are  3-11,  11-15,  and  15  -  15  x  10  ^ 

3  -1  -1 

cm  molec  sec  for  the  same  v  ranges.  The  vibrational  energy  loss  upon  electronic 

state  transfer  was  best  described  by  a  Poisson  type  function  displaced  to  lower 
energy  from  the  initial  energy.  These  basic  models  also  describe  the  relaxation 
in  the  other  gases  with  He  and  Ne  being  less  efficient  and  Kr  and  more 
efficient  than  Ar.  The  magnitudes  of  the  rate  constants  and  the  models  are 


discussed. 
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I.  INTRODUCTION 

Because  of  the  importance  of  rare  gas  halide  excimer  lasers,  considerable 

research  has  been  directed  toward  characterizing  and  understanding  the  basic 

kinetic  processes  of  formation,  quenching,  electronic  state  transfer,  and 

vibrational  relaxation  of  the  B ( 1 / 2 )  and  C ( 3 / 2 )  states  of  the  rare  gas  halide 

molecules.'*’  The  importance  of  the  C  state  has  only  recently  been  realized.  In 

2 

fact,  the  C-A  transition  of  XeF  has  been  made  to  lase.  Clt  irly  the  C  state  must 

be  included  in  any  realistic  model  of  the  laser  medium  and  understanding  the  B-C 

state  coupling  and  the  associated  vibrational  relaxation  may  be  of  importance  in 

optimizing  the  design  and  performance  of  the  rare  gas  halide  lasers. 

Previous  reports  from  this  and  other  laboratories  have  provided  information 

3 

concerning  the  reactions  of  the  first  excited  states,  Rg(  P0  ),  of  the  rare  gas 

atom  with  halogen  donating  molecules  to  form  the  RgX(3,C)  molecules.  Ouenching 

3  3 

rate  constants  and  branching  fractions  have  been  report  d  for  Ar(  Kr  (  P?), 

3  3_4 

and  Xe(  P^)  with  a  variety  of  halogen  donors.  In  addition,  the  bound-free 

emission  spectra  from  the  RgX(B-X)  and  RgX(C-A)  transitions  have  been  interpreted 

to  provide  information  concerning  the  vibrational  energy  disposal  and  dynamical 

features  of  these  reactions. ^ ^  Although  most  of  the  work  has  been  concerned 

•  3  3  7 

with  the  metastable  ?2  state,  results  for  the  P^  state  appear  quite  similar. 

In  this  paper,  we  present  results  concerning  the  vibrational  relaxation 

and  electronic  state  transfer  for  the  XeCl(B,C)  states  in  various  bath  gases. 

We  chose  XeCl  because  of  the  ease  of  obtaining  experimental  data  and  because 

assignment  of  vibrational  distributions  from  the  bound-free  spectra  is  more 

straightforward.  The  general  findings  from  the  XeCl  results,  should  be  extendable 

* 

to  other  RgX  systems.  We  discuss  in  detail  the  results  for  XeCl  in  argon  and  sub¬ 
sequently  contrast  the  relaxation  processes  in  Ar  with  those  in  P.e,  Ne ,  Kr ,  Xe 
and  N^.  The  final  results  are  in  the  form  of  state-to-state  rate  constants  for 
C-B  transfer  and  for  vibrational  relaxation  as  a  function  of  vibrational  level. 

In  addition  to  applications  to  rare  gas  halide  lasers,  the  results  are  of  general 


interest  for  relaxation  of  high  vibrational  levels  of  ionic  diatomic  molecules 
in  rare  gases. 

3 

Reactive  quenching  of  Xe(  P?  with  halogen  donors  gives  XeCl(B,C)  with 
various  vibrational  distributions  and  using  different  donors  permits  a  range  of 
the  vibrational  levels  to  be  emphasized.  For  instance,  the  donors  Cl?,  CCl^  and 
C0C1-,  give  vibrational  distribution  peaking  at  v=0,  50  and  100,  respectively.^ 

Computer  simulation  of  the  B-X  and  C-A  emission  spectra  permits  assignment  of  the  vibra 
tional  distribution.  Initial  distributions  are  obtained  from  low  pressure  (<0.2 
torr)  spectra  where  collisions  with  the  bath  gas  are  net  important.  At  higher  pressures 
the  appearance  of  the  B-X  and  C-A  emissions,  as  well  as  their  relative  intensities, 
change  as  a  consequence  of  vibrational  relaxation  and  transfer  between  the  two 
electronic  states  by  bath  gas  collisions.  By  modeling  the  population  distribu¬ 
tions  and  the  B/C  intensity  ratios  as  a  function  of  bath  gas  pressure,  rate  constants  f 
relaxation  and  transfer  can  be  deduced,  as  well  as  a  state-to-state  view  of  these 
processes.  The  state-to-state  kinetics  are  emphasized  in  this  paper;  however,  the 
results  are  of  practical  relevance  because  the  rare  gas  halide  formation 
processes  in  lasers,  whether  they  be  reactive  quenching  of  excited  rare  gas  atoms 
or  ion-ion  combination,  initially  yield  high  vibrational  levels  of  RgX(B,C) . 

The  organization  of  the  paper  is  as  follows.  In  section  II  the  experimental 

3 

results  for  the  reactions  of  Xe(  with  C^,  CCl^,  and  COCl^  in  various  bath 

gases  are  summarized.  Also  the  XeCl(B,C)  vibrational  distributions  obtained  from 
simulation  of  spectra  are  presented.  The  XeCl(B-X)  and  XeCl(C-A)  simulation  is 
nearly  the  same  as  used  previously^  this  is  only  briefly  mentioned.  Section 

III  gives  the  details  of  the  modeling  procedures  used  to  reproduce  the  variation 

of  the  populations  with  pressure.  The  detailed  results  for  Ar  are  given  in  Section 

IV  and  a  summary  of  the  other  gases  is  given  in  Section  V.  A  discussion  is  given 
in  Section  VI. 


II.  EXPERIMENTAL  RESULTS  AND  SPECTRAL  SIMULATIONS. 
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The  experimental  details  for  stuu.  of  the  Xe(  P?)  and  Xe(  P  )  atom  reaction 
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with  chlorine  donors  have  been  given  elsewhere  ’  and  can  be  briefly  summarized 

here.  The  metastable  atom  reactions  are  studied  with  a  flowing  afterglow  apparatus 

A  small  amount  of  Xe  is  entrained  in  an  Ar  carrier  gas  flow  which  is  passed  through 

3 

a  low  power  hollow  cathode  discharge.  Xe(  P^)  is  formed  from  direct  excitation 

3 

in  the  discharge  and  from  excitation  transfer  from  Ar(  P  )  +  Xe.  Only  the  rela- 

3 

tively  long  lived  Xe(  P^)  survives  the  transit  time  to  the  reaction  zone  several 
cm.  downstream  of  the  discharge.  The  chlorine  donor  is  added  coaxially  to  the 

X 

Ar/Xe  flow  in  the  reaction  zone  and  the  resulting  XeCl  emission  is  observed 
with  a  monochromator  and  detection  system  interfaced  to  a  minicomputer.  The  com¬ 
puter  controls  data  acquisition  and  stores  the  spectra  on  magnetic  tape  for  sub- 

g 

sequent  analysis.  The  spectra  can  be  corrected  for  the  spectral  response  of  the 

detection  system,  integrated  to  obtain  (B-X)  and  (C-A)  emission  intensities,  and 

replotted  for  direct  comparison  with  simulated  spectra.  The  pressure  at  which 

the  spectra  are  observed  can  be  changed  by  varying  the  Ar  flow  and/or  throttling 

the  pumping  speed.  The  flowing  afterglow  method  is  somewhat  limited  in  that  only 

inexpensive  gases  (Ar)  can  be  routinely  used  as  the  carrier  gas  and  the  pressure 

range  is  restricted  from  0.2  to  ^  50  t'orr. 

For  the  resonance  state  atom  reactions  a  mixture  of  Xe,  chlorine  donor,  and 

bath  gas,  under  either  static  or  flowing  conditions  is  irradiated  in  a  reaction 

cell,  with  a  microwave  powered  Xe  resonance  lamp  interfaced  to  the  cell  with  a 

MgF^  window.  R.adiation  trapping  in  the  cell  lengthens  the  effective  lifetime  of 
3 

Xe(  P^)  to  a  time  scale  permitting  reaction  to  occur.  The  XeCl*  emission  is  ob¬ 
served,  at  right  angles  to  the  light  from  the  resonance  lamp,  with  the  same  de¬ 
tection  system  as  above.  Relative  to  the  flowing  afterglow,  this  method  has  the 
advantage  of  a  much  larger  pressure  range,  0.1  torr  to  5  atm,  as  well  as  permitting 
different  bath  gases  to  be  used.  One  disadvantage  is  scattered  light,  mainly  Xe* 
lines,  from  the  discharge  which  overlap  the  long  wavelength  part  of  the 
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XeCl(C-A)  spectrum.  The  intensity  of  the  atomic  lines  must  be  carefully  subtracted  out 
to  obtain  accurate  B/C  intensity  ratios.  In  all  of  the  work  here  the  small  (the 
calculated  B-A/B-X  branching  ratio  is  0.06)  contribution  of  the  B-A  transition  to 
the  C-A  spectrum  has  been  ignored. 

Sample  XeCl  spectra  obtained  with  the  above  methods  are  shown  in  Figs.  1-3 

for  various  donors  and  pressures  of  argon.  Although  close  examination  does  show 

slight  differences  resulting  from  the  978  cm  ^  higher  energy  for  the  resonance 
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state,  the  spectra  from  the  Xe(  P^)  and  Xe(  ?^)  reactions  are  virtually  identical 
and  in  the  remainder  of  this  paper  no  distinction  will  be  made  between  them.  Plots 
of  the  intensity  ratio,  I  /I  ,  vs  pressure  are  shown  in  Fig.  4a-b  for  the  reactions 
of  Xe*  +  Cl-,  CCl^  and  COC^  in  argon  and  for  Xe*  +  Cl-  in  various  bath  gases. 

The  numerical  simulation  of  the  experimental  bound-free  spectra  has  been 
presented  elsewhere.^  Three  basic  variables  enter  into  the  simulations;  relative 
shapes  of  the  two  potentials,  variation  of  the  transition  moment  with  inter- 
nuclear  distance  and  the  vibrational  distributions.  The  variation  of  the 
transition  moment  is  especially  important  in  fitting  the  bound-free  emission.  In 
this  work,  we  used  the  procedure  developed  previously  for  fitting  the  XeCl(B-X) 
emission,  which  employed  Tellinghusin' s  B  state  potential  in  a  Rittner  form  with 
a  flat  potential  for  the  X  state. ^  The  effects  of  the  transition  moment  was 
included  by  introducing  a  monotonic  function  to  multiply  the  spectrum  obtained 
using  a  constant  transition  moment.  These  results  were  successful  in  simulating 
the  XeCl(B-X)  emission  from  a  variety  of  donors.  The  B-X  transition  is  easy  to 
simulate  because  the  spectra  are  not  very  sensitive  to  the  shape  of  the  lower 
potential;  such  is  not  the  case  for  the  C-A  transition  where  the  repulsive  wall 
for  the  lower  state  A-potential  is  very  important.  Initially  we  used  the 
potentials  previously  employed  to  simulate  the  C-A  transition.  However,  when 
modeling  the  XeCl(B)  distribution  as  a  function  of  pressure,  a  second  maximum 
developed  in  the  distribution  from  C-*B  collisional  transfer.  This  second  peak, 
which  is  not  experimentally  observed,  was  a  consequence  of  the  initial  C  distribution 
having  its  maximum  at  lower  E,,  than  for  XeCl (3) .  This  is  evident  in  Fits.  -  and  7 
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of  ref  5b  where  E  (XeCl.B)  is  6  kcal  mole  ^  larger  than  E  (XeCl.C)  from  the 
—  mp  mp 

Xe*  +  Cl?  reaction.  To  rectify  the  problem  it  was  necessary  to  make  the  A 
state  slightly  more  repulsive  so  that  the  vibrational  distribution  from  the  C-A 


emission  was  moved  to  higher  E^.  The  initial  XeCl(C)  distributions  used  here 

from  Cl 2 ,  CCl^  and  C0C19  were  indistinguishable  from  the  XeCl(B)  distributions. 

The  expressions  for  the  potentials  used  here  to  simulate  the  C-A  emission  are 

C ,  V  2 


V(R) 


C  »3/2  He.  1 - vT(R') 

B,1/2V  U  ' 


B.l/2 


cm 


V(R)A,3/2  =  1>7?  x  1Q7  exp(_2.827  R') 


cm 


a) 


b) 


(1) 


where  V^(R')R  is  given  by  eqn  4a  of  ref.  5b_,  R'  =  +  0 . 9814  (R-R  ^  ’ 3^2)  - 

0.0365(R-ReC,B/2)2,  ^eB,1/2  =  195.2  cm-1,  u^0,372  =  188  cm"1,  ReB’1/2  =  2.937A, 

C  3/2  0 

and  R  ’  =  3.14A.  It  should  be  noted  that  the  expressions  for  the  A  potential 


and  the  transition  moment  function  of  ref.  5b  (eqn  14  and  15)  are  slightly  in 
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error. 

The  basic  procedure  for  simulating  the  experimental  spectra  is  to  first 
calculate  spectra  for  individual  vibrational  levels  of  the  upper  state.  These 
individual  spectra  are  then  summed  over  an  assumed  vibrational  distribution  to 
obtain  a  total  spectrum.  Different  vibrational  distributions  were  tried  until 
a  satisfactory  fit  was  found  for  the  experimental  spectra  obtained  at  different 
pressures.  Some  simulated  spectra  are  compared  to  the  experimental  ones  in 
Figs.  1-3  and  the  corresponding  steady-state  vibrational  distributions  are  shown 
in  Fig.  5-7. 


Two  important  points  affecting  the  modeling  are  the  reliability  of  the 
distributions,  i.e.  how  accurate  are  the  simulations,  and  the  variation  of  the 
XeCl*  lifetime  with  v  level.  The  first  point  is  addressed  in  Fig.  8,  which  shows 
XeCl*(B-X) emission  spectra  in  Ar  from  Xe*  +  Cl0  at  1.3  torr  and  from  Xe*  +  CCl^ 
at  11.9  torr  along  with  two  simulated  spectra  that  were  judged  to  be  satisfactory 
fits.  The  corresponding  vibrational  distributions  differ 


by  ^  10%  on  Che  average.  Hence,  any  relaxation  model  that  matches  a  distribution 

to  within  'v  10%  of  the  distribution  assigned  by  simulation  must  be  considered 

satisfactory.  The  uncertainty  in  the  distributions  increases  with  increasing 

pressure  (reduction  in  vibrational  energy)  because  of  the  loss  of  the  oscillations 

from  the  spectra.  This  is  evidenced  by  the  spectrum  from  CCl^  for  which  the 

distributions  are  less  unique  than  for  the  low  pressure  Cl^  spectra.  Above 

■v  15  torr  of  Ar,  the  oscillations  have  disappeared  and  the  best  that  can  be 

done  is  to  assign  high  temperature  Boltzmann  distributions.  The  above 

discussion  was  based  on  the  B-X  spectrum.  For  the  C-A  spectra,  the  potential 

curves  are  less  accurate  and  the  B-A  and  C-A  emission  overlap;  therefore,  the 

distributions  from  the  spectral  simulations  are  less  reliable  than  those  from 

the  B-X  simulations.  Hence  in  modeling,  emphases  was  placed  on  fitting  first 

the  I  /I  ratio  vs  pressure,  the  B  state  vibrational  distribution,  and  finallv 
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the  C  state  vibrational  distribution. 

The  second  point  was  the  variation  of  the  XeCl(B,C)  lifetimes  with  v  level. 

Unfortunately  measurements  of  the  XeCl(B,C)  lifetimes  have  not  been  reported. 

10a 

However,  the  theoretical  values  given  by  Dunning  and  Hay  for  XeCl(B,v=0), 

11  nsec,  and  XeCl(C,v=0),  120  nsec,  can  be  combined  with  the  dependence  of  the 

transition  moment  on  internuclear  distance  to  obtain  calculated  lifetimes. 

The  Einstein  coefficients  decrease  by  ^  3  from  v=0  to  100;  therefore,  the  lifetimes 

increase  by  'v  3.  This  trend  of  increasing  lifetime  vs.  v  level  has  been  observed 

for  a  small  range  of  v  for  KrCl. ^  Dunning  and  Hay  have  calculated  lifetimes  for 

ArF,  KrF  and  all  the  xenon  halides. ^  Experimental  lifetimes  have  been  reported 

for  KrFCB)^3  and  XeF(B,C)^k;  the  results  tend  to  agree  with  the  calculated 

values  to  within  'v  30%.  The  XeCl  calculated  lifetimes  and  energy  levels  used 

in  the  modeling  calculations  are  summarized  in  Table  1. 

Even  without  the  modeling  results,  some  general  conclusions  can  be  drawn  from 
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the  data.  ’  The  statements  below  refer  to  argon  bath  gas  but  the  same  trends  are 


observed  in  other  bath  gases;  but,  the  characteristic  pressure  regions  are 

different.  At  zero  pressure,  the  initial  I  /I  is  1.3  to  1.5  for  all  3  donors, 

Cl?,  CCl^,  COC^.  In  the  low  pressure  (<_  5  torr)  regime,  B  state  vibrational 

relaxation,  C  state  vibrational  relaxation,  and  a  rapidly  increasing  Ig/I  is 

observed.  Since  at  these  pressures,  the  lifetime  of  the  B  state  is  less  than 

the  time  between  collisions  with  bath  gas,  the  apparent  B  state  relaxation  is 

due  to  collisional  transfer  from  the  C  state.  This  is  also  the  cause  of  the 

increasing  I  /I  ratio.  A  simple  steady  state  analysis  for  XeCITB)  and  XeCl(C)  similar 

A  3.  7 

to  those  used  previously  ’  can  be  used  to  obtain  approximate  rate  constants  for  C 
to  B  transfer,  k  ,  from  the  variations  of  I  /I  vs  pressure  in  the  low  pressure  regime 
The  analysis,  which  included  an  average  radiative  lifetime  estimate  for  each 
electronic  state  gave  kcg  of  18,  17,  5.9,  and  2.7  x  10_11  cm3  molec-1  sec-1 
for  Kr,  Ar,  Ne,  and  He, respectively .  Since  the  reaction  of  Cl^  with  Xe*  gives 
a  high  vibrational  distribution,  these  rate  constants  refer  to  v  v  100.  The  Ig/I 
ratios  from  CC1,  and  C0C1„  do  not  increase  as  fast  as  for  Cl„ ;  hence,  the  transfer 
rate  constants  evidently  are  smaller  for  lower  v  levels.  This  is  supported  by  the 
fact  that  the  change  in  Ig/Ic  "slows  down"  for  higher  pressure  and  in  the  inter¬ 
mediate  regime  (5  to  50  torr  of  Ar)  I  /I  is  nearly  constant.  In  this  regime  coll- 
isions  with  the  B  state  become  important  and  the  slow  change  of  I  /I_  is  due  a 
combination  of  reduced  transfer  rate  constants  for  the  lower  v  levels  and  also  to 
the  balance  in  B  and  C  state  populations  as  the  B-C  transfer  rate  competes  with  radiative 
decay  of  the  B  state.  At  v5C  torr  the  I  /I  attains  its  maximum  value.  Examination 

of  the  10-100  torr  spectra  indicates  that  the  vibrational  distributions  are  high 
temperature  Boltzmann  distributions,  while  the  electronic  state  distribution  is 
non- Boltzmann.  By  v  100  torr,  the  vibrational  distributions  are  nearly  300K.  As 
the  pressure  is  increased  to  the  1000  Torr  region,  the  I  /I  ratio  decreases  as 
transfer  between  B  and  C  occurs  to  give  a  Boltzmann  electronic  state  distribution. 

At  the  highest  pressures  electronic  transfer  is  complete  and  interpreting  the  high 


pressure  I  /I  ratio  as  a  300K  distribution  shows  that  the  XeCl(C)  state  is 
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'v.  200  cm  ^  lower  in  energy  than  the  XeCl(B)  state.'7  The  exact  value  of  the  energy 

1  ^3 

separation  is  still  being  discussed.  Tellinghusien  and  McKeever  report  a  value 

of  'v  130  cm  ^  by  interpreting  the  X  / 1^  high  pressure  intercept  at  360K  and 

B  C  A 

c _A 

800K.  (In  addition  they  were  able  to  obtain  - -  'v  t  /t  of  0.087  which  is 

^-X  14b 

excellent  agreement  with  the  calculated  value  of  0.091).  Bokor  and  Rhodes  , 

on  the  other  hand,  report  an  energy  separation  of  ^  5+25  cm  ^  with  the  state 

higher  in  energy.  However,  their  experiments,  which  used  an  ArF  laser  to  excite 

Xe/C^  mixtures  in  Ar,  extend  to  only  1000  torr.  In  fact  pressures  of  >  1000 

torr  are  required  to  obtain  B/C  equilibrium  concentrations  (see  Fig.  2  of 

ref.  7  for  example).  Thus,  the  data  of  Bokor  and  Rhodes  do  not  extend  to 

sufficiently  high  enough  pressure  to  obtain  a  true  equilibrium  I_/I_  ratio. 

D  C 

Experiments  of  the  type  described  in  ref.  7  have  recently  been  performed  at  low 
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temperatures  and  preliminary  results  support  an  energy  separation  of  ^  200 

cm  .  The  uncertainty  of  the  B-C  energy  separation  does  not  affect  interpretations 

of  this  paper  because  we  are  concerned  with  collisions  of  XeCl  molecules 

in  high  v  levels. 


III.  METHODS  FOR  MODELING  VIBRATIONAL  RELAXATION  AND  B-C  STATE  TRANSFER 

The  method  used  to  model  the  variation  of  the  XeCl(B,C)  vibrational  distri¬ 
butions  and  B/C  ratio  as  a  function  of  pressure  is  essentially  the  same  as  the 


one  used  by  Duewer,  et  alx  in  studying  the  rotational  relaxation  of  CN(B,v=ll). 
The  one  important  difference  here  is  the  large  variation  in  radiative  lifetime 
which  must  be  included  in  the  calculation.  Electronic  quenching 


of  XeCl(C,B)  is  assumed  to  be  unimportant  since  the  modeling  only  extends  to 
50  torr.  The  vibrational  levels  of  XeCl(B,C)  are  divided  into  three  ranges; 
low  v  (v=0-30) ,  mid  v  (v=30-70),  and  high  v  (v=70-130) .  Examination 


of  Figs.  5-7  shows  that  Xe*  +  gives  an  initial  vibrational  distribution 
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peaked  at  v  ^  100  with  _>  80%  of  the  total  in  the  high  v  range;  Xe*  +  CCl^  gives 
an  initial  distribution  peaked  at  v  'v  50  with  80%  in  the  mid  v  range,  and  COCl0 
gives  a  distribution  marked  at  v=0  wii_h  80%  in  the  low  v  range. 

If  the  relative  populations  of  XeCl(B,v)  and  XeCl(C,v)  are  considered  to 
be  a  vector,  then  the  normalized  initial  population  vector,  X(0),  is  given  by 

X(0)  =  R  /ER  (2) 

where  R^  is  the  formation  rate  into  the  ith  level.  If  NC  and  NB  are  the  maximum 
number  of  vibrational  levels  for  XeCl(C)  and  XeCl(B),  respectively,  the  total 
C  and  B  populations  are 

NC  '  NC+NB 

N~  *  E  X.(0)  and  Na  =  £  X.(0).  (3) 

C  i=l  1  B  NC  1 

Since  the  reaction  of  Xe*  +  C^  populates  levels  up  to  v=140  in  the  C  state, 

X(0)  would  include  over  200  levels  if  each  vibrational  level  was  included. 

To  save  computational  time,  the  populations  were  grouped  into  every  4th  vibrational 
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level.  The  relative  B  to  C  formation  rate  is  given  by  the  zero  pressure  I  /I 

B  C, 

ratio,  which  was  taken  as  1.4  for  all  three  donors  from  Fig.  4.  The  initial 
relative  vibrational  distribution  within  each  electronic  state  was  obtained  from 
the  simulation  of  the  lowest  pressure,  typically  £  0.2  torr,  spectrum.  The  low 
pressure  distributions  shown  in  Fig.  5-7  are  steady-state  distributions  and 
must  be  corrected  for  the  lifetime  variation  to  obtain  the  true  initial  distri¬ 
butions.^ 

The  relative  populations  that  exist  after  exactly  m  collisions  are  given  by 

X(m)  =  ?m  X(0)  (4) 

The  transition  matrix  has  elements,  P ,  which  give  the  probability 
per  collision  with  bath  gas  of  a  XeCl  molecule  in  state  j  going  to  state  i.  P  was 
partitioned  into  4  parts  as  follows: 


for  1  <.  i  <_  NC 
for  1  <  i  <  NC 
for  NC+1  ^  i  +  NC+NB 
for  NC+1  <  i  <  NC+NB 


1  <  j  <  NC 
1  <  j  <  NC  +  NB 
1  <  j  <  NC 
NC+1  <  j  <  NC+NB 


p?. 

=  c 

state 

relaxation 

pt. 

ij 

=  B 

to  C 

transfer 

pT. 

ij 

=  C 

to  B 

transfer 

pR. 

ij 

=  B 

state 

relaxation 

The  diagonal  elements  of  the  matrix  correspond  to  elastic  collisions,  which  are 
unobservable  in  this  case,  and  were  calculated  from  the  normalization  condition. 


P..  =  1-IP..  i^j  (6) 

J 1  ±  iJ 

By  normalizing  ^  in  this  manner,  the  sum  of  the  inelastic  and  elastic  rate 
processes  is  independent  of  vibrational  level  providing  the  collision 
frequency  (cross-section)  is  selected  to  be  constant  for  all  levels.  In  the 
modeling  the  collision  frequency  was  taken  as  a  constant;  but,  the  ratio  of  the 
inelastic  and  elastic  cross-sections  was  allowed  to  be  a  function  of  vibrational 
level.  The  transition  matrix  is  subject  to  detailed  balance  and  up  transitions  were 
calculated  from  the  model  assumed  for  the  down  transitions,  i.e.  for  j  >  i. 

•  P.,  ,-W*  (7) 


P..  .  P..  e-(ErEl)/ltT 

jl  1J  J 


The  probability  that  a  molecule  undergoes  exactly  m  collisions  before  decaying 
by  photon  emission  can  be  calculated  as  follows.  For  a  molecule  in  a  given  level 
with  a  radiative  lifetime,  x,  the  distribution  of  times  for  radiative  decay  is 
given  by 

F(t)  =  t-1  e~t/T  (8) 

The  fraction  of  molecules  undergoing  exactly  m  collisions  in  a  fixed  time 
interval  t  is  given  by 

F(m)  =  e"Zt  (9) 

m ! 

where  z  is  the  collision  rate  corresponding  to  a  given  pressure.  The  distribution 
for  the  number  of  collisions  before  radiative  decay  is  given  by  the  convolution 


of  these  two  distributions  over  all  time, 
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where  A(m)  is  a  vector  and  is  different  for  each  level  because  of  the  different  values 

of  Although  this  equation  contains  an  infinite  number  of  terms,  the  summation 

can  be  approximated  with  a  finite  number  of  terms,  m*,  because  A(m)  decreases 

with  increasing  m.  In  this  work  the  summation  was  terminated  when  EX.(m*)  was 

i 

0.1%  of  X^(0).  For  this  criteria  m*  usually  ranged  from  2  to  100  depending 
on  pressure  and  bath  gas.  The  intensity  from  individual  levels  is  summed  to 
obtain  the  Ig/Ic  ratio.  The  intensities  were  converted  to  steady-state  concen¬ 
trations,  X^(s) ,  by 

Xi(s)  -  I.t.  •  (12) 

If  lifetimes  for  all  levels  were  equal,  the  relative  steady-state  concentrations 
and  intensities  would  be  equivalent  and  no  conversion  would  be  needed.  The 
distributions  given  by  X^Cs)  are  compared  to  the  distributions  obtained  from 
simulation  of  the  experimental  spectra  at  a  given  pressure  to  gain  information 
about  the  form  of  P .  . . 

ij 

Assignment  of  the  model  for  £  was  based  on  a  trial  and  error  approach.  We 
began  with  simple  models  and  added  more  flexibility  as  necessary.  We  wanted 
overall  rate  constants  for  transfer  and  vibrational  relaxation,  as  well  as  state- 
to-state  rate  constants.  The  state-to-state  processes  are  determined  by  the 
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form  of  ?  and  the  rate  constants  by  zip  (i^j).  The  success  of  a  model  for  P  is 

i  ^ 

judged  by  how  well  the  calculated  distributions  fit  the  variation  of  the  shape 
of  the  simulated  distributions  with  change  in  pressure.  On  the  other  hand,  the 
rate  of  relaxation  of  the  vibrational  distributions  as  a  function  of  pressure 
depends  on  both  the  collision  rate  z  (or  a,  the  total  collision  cross-section)  and 
For  our  purposes  the  only  restriction  for  o  is  that  it  be  large  enough  to 
give  the  required  rates  of  relaxation  and  transfer.  Initially,  we  estimated 
"reasonable"  values  for  the  collision  cross-section  and  proceeded  with  the  modeling. 
But  in  some  cases  it  was  subsequently  necessary  to  increase  a  because  the  required 
rates  of  relaxation  and  transfer  could  not  be  made  fast  enough  for  the  first  choice 
for  a.  After  selection  of  a  P  that  fitted  the  overall  shape  of  the  distributions, 
though  not  necessarily  at  the  right  pressures,  the  inelastic  cross-section  could  be 
increased  or  decreased  until  the  distributions  matched  at  the  proper  pressure. 

The  objective  is  to  obtain  a  model  for  the  transition  matrix,  P,  which 
simultaneously  fits  all  of  the  experimental  data;  i.e.  the  B  and  C  distributions 
and  the  Ig/I^,  ratio  as  a  function  of  pressure  for  all  three  Cl  donors.  Preferably 
the  model  should  be  described  by  as  few  independent  parameters  as  possible.  Cal¬ 
culations  first  were  done  for  XeCl  in  Ar.  Relaxation  and  transfer  are  both 
sufficiently  fast  that  changes  occur  at  relatively  low  pressure.  Since  the  required 
number  of  interations  needed  for  convergence  depends  on  pressure,  analysis  of  the 
argon  data  for  a  given  degree  of  relaxation,  requires  less  computer  time  than  for 
He  or  Ne.  The  fitting  for  Ar  also  provides  a  good  reference  point  for  both  lighter 
and  heavier  bath  gases.  To  aid  in  the  discussion  and  in  identification  of  results, 
each  model  will  be  given  a  label  (capital  letter)  as  it  is  presented.  Several 
different  models  were  tried.  Each  of  these  models  has  several  parameters  which 
can  be  varied,  but  we  present  only  the  best  set  of  parameters  from  the  several 
trial  and  error  attempts  to  obtain  a  satisfactory  fit  to  the  experimental  data. 

Table  2  contains  a  summary  of  the  models  and  the  collision  cross-sections. 


IV.  RELAXATION  IN  ARGON 


IV  A.  Vibrational  Relaxation  of  XeCl(C,  high  v) ,  Model  A. 

The  XeCl(B)  radiative  lifetime  is  less  than  the  time  between  collisions  for 
pressures  <_  2  torr  and,  B  to  C  transfer  can  not  be  important.  Therefore,  the 
change  in  the  XeCl(C)  distribution  is  solely  due  to  vibrational  relaxation. 

Although  this  statement  is  qualitatively  correct,  20%  of  the  XeCl(B,v=100)  molecules 


undergo  at  least  one  collision  at  0.5  torr.  This  percentage  increases  to  40%  at  1 
torr  and  'v  55%  at  2  torr;  but,as  a  first  approximation^  to  C  transfer  can  be 
ignored  and  C  state  relaxation  can  be  modeled  as  vibrational  relaxation  for  pressures 
£  2  torr. 

£ 

The  initial  model  used  to  calculate  the  transition  matrix  elements,  P..,  for 
transitions  j  >  i  was  given  by 

pij  =  fr  +  exp(-  (13) 

with  AE  =  E_.  -  E^.  The  constants  a  and  b  were  treated  as  adjustable  parameters 
and  were  varied  to  obtain  a  satisfactory  fit  to  the  data.  The  kT  factor  (k  = 
Boltzmann  constant  and  T  *  300°K)  serves  to  make  AE  unitless  and  has  no  other 


physical  significance.  For  some  values  of  a  and  b,  z  pK  was 

r 

greater  than  1  and  values  of  P^  calculated  by  the  above  equation  were  multiplied 

£ 

by  a  normalization  constant  so  that  E  P..  =  1.  Certain  variations  of  the  above 

i 

model  also  were  explored;  P^.  =  a  +  exp(-bAE/KT) ,  P^  =  akT/AE,  and  P^  = 
exp(-bAE/KT) ,  and  close  fits  were  obtained  for  most  of  these  variations.  However, 


the  "best"  overall  fit  was  obtained  with  eqn.  13  for  a  =  1.0  and  b  =  0.11.  With 
normalization  the  best  model  for  ?  was 

D  1  k  T 

Pij  =  l0  (AE  +  exP<-0-U  AE/kT)),  (14) 

O 

which  will  be  referred  to  as  Model  A,  with  a  collision  diameter  of  4.95  A.  The 

£ 

dominant  contribution  to  the  magnitude  of  P, .  is  given  by  the  exponential  term  and 


since  the  vibrational  relaxation  could  be  fit  satisfactorily  with  just  an  exponential 
term,  the  other  models  presented  later  do  not  include  the  kT/AE  term. 

The  population  distributions  resulting  from  model  A  are  compared  in  Fig.  9 
to  the  simulated  distributions  at  several  pressures.  The  fit  is  well  within  the 
10%  uncertainty  of  the  experimental  distributions  throughout  the  0.11  to  2.2  torr 
range.  Note  especially  the  good  fit  for  P  =  2.2  torr  for  20  <_  v  £  110.  By 


combining  the  collision  frequency  and  the  Z  P  values  the  relaxation  rate  con- 
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stants,  k^(10  cm  molec  sec  )  as  a  function  of  v  level  are  4.5  -  13, 


13  -  20,  and  20  -  34  for  v  =  4-30,  30-70,  and  70-130,  respectively.  The 


smooth  increase  of  the  rate  constants  with  v  is  due  to  the  smaller  energy 


separation  between  high  v  levels  and  the  larger  number  of  levels  below  a  given 


v  level. 

£ 

For  model  A  the  P„  decline  smoothly  with  increasing  AE.  Single  quantum  jumps, 
(which  for  our  mod',-1  includes  all  transitions  from  v  to  v-1,  v-2  ,  v-3,  and  v-4)  are 


the  most  important;  however,  multiquantum  transitions  also  are  significant.  The 
adjustable  parameter  b  in  the  exponential  term  governs  the  range  of  AE.  For 
small  b,  the  low  v  range  gains  in  population  at  the  expense  of  the  high  v  and 
the  distribution  is  broadened,  whereas  for. larger  b  the  distribution  remains, 
sharp  and  moves  systematically  to  lower  v.  For  a  successful  fit  b  must  be 
0.1  -  0.2.  The  success  of  this  model  is  due  to  the  fact  that  molecules  in  the 
high  v  range  (v  =  70  -  130)  have  a  significant  probability  of  transfer  to  the 
low  v  range  (v  *  4  -  30).  This  must  be  a  general  feature  of  all  models  that 
fit  the  relaxation  of  XeCl(high  v)  in  Ar. 

Since  the  distribution  from  Xe*  +  Cl^  is  in  the  high  v  range,  the  model 
should  be  most  reliable  for  these  levels.  Most  of  the  population  in  the  low  v 
levels  from  Xe*  +  Cl^  arises  from  multiple  collisions  and  the  model  may  not  be 


too  reliable  for  the  transition  probabilities  of  the  lower  levels. 


IV-B.  Vibrational  Relaxation  and  Transfer  of  XeCl  (high  v)  -  Model  B 

Before  proceeding  with  the  addition  of  the  XeCl(B)  state  to  the  calculations, 
certain  assumptions  must  be  made.  Because  C  to  B  transfer  is  an  important  process 
at  low  pressure,  apparent  vibrational  relaxation  of  B  is  observed  before 
bath  gas  collisions  with  B  are  important.  At  higher  pressures  where  XeCl(B)  + 
bath  gas  collisions  become  important,  the  vibrational  distribution  of  XeCl(B) 
has  relaxed  so  much  that  obtaining  information  about  XeCl(B)  vibrational  re¬ 
laxation  is  impossible.  Without  any  knowledge  to  the  contrary,  we  will  assume 
that  B  state  relaxation  can  be  described  by  the  same  model  as  for  XeCl(C). 

We  will  drop  the  kT/AE  term  because  of  its  minor  role,  and  the  model  for  vibrational 
relaxation  becomes 

pj.  =  exp(-bAE/kT) .  (15) 

Since  the  vibrational  levels  of (C)  are  closer  together  than  for  (B),the  k  (B)  are 
slightly  smaller,  typically  by  10-20%,  than  k^CC)  for  the  same  v  level. 

The  additional  assumption  concerns  B  ->•  C  transfer.  For  the  same  reasons 
as  given  above,  reliable  information  for  B-C  transfer  from  XeCl(B,  high 
v)  is  not  available.  However,  detailed  balance  can  be  used  to  relate  the  forward 
and  reverse  B-C  transfer  rates.  At  limiting  high  temperature  the  populations  of 

both  states  wouia  be  equal  and  k  =  k__.  Although  the  translational  temperature 

CB 

is  near  300K,  it  is  reasonable  to  assume  that  for  high  vibrational  excitation 

kD_(v)  v  k__(v)  and  that  the  form  for  P^.  is  the  same  for  both  XeCl(B)  and  XeCl(C). 
ot  CB  ij 

The  XeCl(C,v=0)  level,  which  represents  the  total  v=0-3  population)  is  the  lowest 
level  and  has  no  down  transition,  hence,  its  kT  and  P  values  are  determined  completely 
from  detailed  balance.  The  course  graining  will  result  in  an  increasing  erroneous 


representation  as  the  XeCl  population  builds  up  in  the  lowest  levels  of  B  and  C. 
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In  contrast  to  modeling  the  C  state  vibrational  relaxation,  simultaneous 
fitting  of  the  B  and  C  distributions  and  the  I  /I  ratio  as  a  function  of  pressure 

D  L 

was  extremely  difficult.  Considerable  effort  was  expended  in  obtaining  a  model 

T 

that  was  sufficiently  flexible.  We  initially  used  the  same  form  for  P..  that 

ij 

was  used  for  vibrational  relaxation  and  varied  the  parameters,  as  well  as  the 
relative  importance  of  transfer  vs.  relaxation.  For  a  given  C  level  the  highest 
transfer  probability  will  be  to  the  adjacent  lower  B  level,  according  to  this 
model.  For  certain  values  of  the  parameters,  the  I  /I  ratio  could  be  matched 

D  C 

(and  also  C  state  vibrational  relaxation)  ,  but  the  B  state  vibrational  distri¬ 
bution  was  always  too  high  in  energy.  An  additional  problem  arises  because  the 
energy  difference  between  adjacent  B  and  C  levels  (note,  our  calculation  uses 

every  fourth  real  level)  is  not  constant  with  v,  because  the  values  of  m  and 

e 

T 

0>exe  for  XeCl(B)  and  XeCl(C)  differ  and  the  calculation  gave  PT.  of  widely 
varying  magnitude  for  adjacent  levels,  which  led  to  spikes  in  the  B  distribution 


which  were  difficult  to  compare  with  the  smooth  simulated  distributions.  Improve¬ 


ments  were  necessary  to  rectify  these  problems. 

T 

The  final  model  adopted  for  P  requires  that  the  most  probable  transition 

from  a  given  J  be  to  the  energy  region  Eq  and  has  the  following  form 

T  kT  E  , 

P.,  =  —  f  i+1  exp(-C-E*/kT)dE*  E*  =  i AE-E  (16) 

i-J  o  o 

where  kT  and  AE  are  defined  as  before  and  C  is  the  adjustable  parameter  and  E  is 

o 

given  by 

E  =  (1.0  -  F*  E'.)  E.  cm-1  (17) 

where  E^l  =  Ej/16000  cm  ^  (16000  cm  ^  is  approximately  the  energy  of  the  highest 
populated  levels  in  the  Xe  +  reaction)  and  f  is  an  adjustable  parameter. 

Eqn.  16  allows  the  most  probable  transfer  to  be  fixed  at  Eq,  with  Eq  being  some 
fraction  of  initial  energy.  With  the  benefit  of  hindsight  the  E^/16000  cm-1 
term  was  added  so  that  the  lower  v  levels  would  lose  less  energy  than  higher  v 
levels.  The  use  of  the  integral  in  eq .  16  eliminates  the  computational  artifact 


of  widely  varying  P_  for  adjacent  levels  and  gives  smooth  steady-state  B  distri¬ 


butions.  However  the  integral  also  gave  nearly  constant  1  P..  values,  which  are 

i  1'1 


independent  of  j,  and  hence,  the  are  independent  of  vibrational  level.  To 


regain  the  ability  to  control  the  magnitudes  of  k^  for  different  v  levels  the 


1  P..  values  were  fixed  bv  a  series  expansion 
i  U 


:  p  =  tl(t,  +  t3  e:  +  t4e:~) 

i<j  ■J  “  J  3 


(18) 


where  T  are  adjustable  parameters.  To  increase  flexibility,  a  similar 


expression  was  used  for  IP.. 

13 


I  PR.  =  R. (R.  +  R,E!  +  R. E ! 2) 

i<j  13  12  3  J  4  J 


(19) 


with  R.  adjustable.  The  factors  R1  and  T1  allow  the  magnitude  of  k  and  k 
1  J.  L  I  R 

to  be  altered  by  changing  the  ratio  of  elastic  and  inelastic  components  of 
the  cross  section.  The  parameters  and  ^2_3'4  Provide  the  functional 


dependence  of  I  P..  on  v  level.  The  most  extensive  experimental  data  exists  for 


Ar  as  the  bath  gas.  For  this  reason  anc*  ^2~^4  were  established  for  Ar 

and  were  not  changed  for  the  other  bath  gases.  Several  calculations  were  done 
to  determine  a  good  set  of  parameters.  Model  B  was  chosen  as  giving  the  best 


fit  to  the  Xe*  +  data.  The  parameters  for  Model  B  are  b  =  c  =  0.1,  f  =  0.15 


and  Ti  =  1.2,  0.01,  0.7,  -0.5,  =  0.1,  0.01,  1.0,  -0.3  for  i  =  1,  2,  3,  4, 


respectively.  The  calculated  results  are  compared  to  the  experimental  distri¬ 


butions  and  I  /I  ratios  in  Fig.  10-11.  The  fit  is  fairly  good  for  the  B-state 
B  C 


distribution.  The  fit  to  the  C-state  distributions,  while  not  as  good  as  for 


Model  A,  is  still  within  the  uncertainty  of  the  distributions.  The  I^/l.  fit 

U  C 


is  excellent  below  5  torr  and  is  only  slightly  low  for  the  higher  pressures. 


In  fact,  the  calculated  I  /I  ratio  begins  to  decrease  at  8  torr  while  the 

B  C 


experimental  points  continue  to  rise.  This  suggests  that  the  transfer  rate  constants 
for  the  lower  v  levels  are  too  low  relative  to  the  high  v  levels  and/or  that  the 


C-state  is  vibrationally  relaxing  too  fast.  In  either  case,  as  the  C-state 


population  builds  up  in  low  levels,  the  C  to  B  transfer  slows  down  with 
increasing  pressure.  The  second  reason  is  substantiated  by  the  8.0  torr  results 
of  Fig.  10,  which  shows  that  the  calculated  populations  are  too  high  in  the  low 
v  range  where  the  C-B  transfer  rate  constants  are  the  smallest.  Since  the 
calculations  match  lower  pressure  data  where  the  distribution  is  largely 
above  v=7 0,  this  model  was  considered  successful  in  fitting  the  Xe*  +  Cl^  data. 

In  the  next  section  Model  B  is  modified  slightly  to  give  a  good  fit  to  the 

Xe*  +  CCl^  data  with  a  distribution  that  initially  is  in  the  mid  v  range  and  the 
question  concerning  the  relaxation  and  transfer  rates  for  v=30-70  will  be 
addressed  then. 

Combination  of  the  transition  probability  matrix  elements  with  the  collision 

frequency  gives  the  rate  constants  listed  in  Table  2.  The  relaxation  rate 

constants  agree  favorably  with  the  results  from  Model  A.  The  k^  calculated  here 

are  slightly  larger  than  for  Model  A  and  as  discussed  above,  may  be  too  large 

for  the  low  and  mid  v  range.  Nevertheless,  the  good  fit  to  the  data  suggest 

that  the  kR  for  the  high  v  range  are  20-35  10  ^  cm“*  molec  ^  sec  ^ .  The  good 

fit  to  the  B  state  distribution  and  the  I  /I  ratio  below  5  torr  suggests 

that  the  k^  are  also  reliable,  at  least  for  high  v  levels.  The  form  of  P 

was  discussed  when  Model  A  was  presented  and  the  results  for  Model  B  are  very 

T  R 

similar  (see  Figure  13);  however,  the  fora  of  P..  differs  from  P, ..  The  C 

2-3  i3 

molecules  enter  the  B  state  in  an  energy  range  that  is  peaked  at  'u  10-15%  below 
the  initial  energy  level. 

IV-C.  Relaxation  and  Transfer  of  XeCl(mid  v)  and  XeCl(low  v)  -  Model  C 

Since  a  satisfactory  fit  was  provided  to  the  Xe*  4-  Cl^  data,  model  B  was 

used  to  examine  the  data  from  the  Xe*  +  CCl^  reaction,  which  emphasizes  the  mid- 

v  range.  The  results  are  shown  in  Fig.  12  for  three  pressures  and,  as  can  be 

seen,  the  B  and  C  states  both  relax  too  fast  and  I  / I  increases  too  slowlv 

B  C 
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with  increasing  pressure  (Fig.  11).  The  following  improvements  are  needed. 

Since  XeCl(C)  relaxes  too  fast,  k^(v=30-70)  should  be  reduced.  The  low  I  /I 

could  be  a  consequence  of  the  (too)  fast  relaxation  of  the  C  state  and/or  k 

being  too  small  for  v=30-70.  Since  the  B  state  also  relaxes  too  quickly,  k 

probably  is  too  large  or  transfer  occurs  with  too  large  a  loss  of  energy.  The 

most  probable  energy  loss  by  transfer  of  molecules  from  C(v=50)  to  B(v)  is 

only  7-8%  of  the  initial  energy  and  further  reduction  of  this  will  not  give  much 

improvement.  Hence,  the  best  way  to  simultaneously  reduce  the  rate  of  B-state 

relaxation  and  increase  the  I  /I  ratio  is  to  increase  the  c  parameter  in  eqn.  16; 

this  gives  a  sharper  v-level  distribution  entering  the  B  state.  In  addition 

decreasing  k  to  fit  the  C  relaxation  will  increase  I  /I  .  The  following 
R  d  L 

modifications  to  Model  B  were  made:  R  was  changed  to  0.50  which  reduced  k 

1  R 

by  a  factor  of  v  2;  the  parameter  C  in  eqn.  16  was  changed  to  0.2  which  reduced 

the  B  state  relaxation  rate,  and  b  was  increased  to  0.15,  which  gave  a  slightly 

better  fit  than  using  b=0.10.  The  rest  of  the  parameters  were  the  same  as  for 

Model  B.  The  results  from  this  calculation,  hereafter  denoted  as  Model  C,  are 

compared  in  Figs.  11  and  12  and  the  fit  is  seen  to  be  much  improved.  The  calculated 

transfer  rate  constants  for  Model  C  are  the  same  as  for  Model  B  and  the  only 

T 

difference  is  the  form  of  P...  The  relaxation  rate  constants  for  Model  C  are 

ij 

1. 0-6.0  and  6.0-12  10  ^  cm^  raolec  ^  sec  ^  for  the  low  v  and  mid  v  range, 

respectively.  The  results  for  low  v  are  less  reliable  than  for  mid  v  because 
low  v  levels  are  populated  by  multiple  collisions.  The  Model  C  rate  constants, 
which  are  approximately  half  the  values  from  Model  B,  should  be  more  reliable 
because  i)  the  fit  to  the  CCl^  results  are  better  and  ii)  the  effect  of  multiple 
collisions  is  less  severe  in  the  mid  and  low  v  range  for  the  CCl^  data  than  for 
the  CI2  data. 

Model  C  was  used  to  examine  the  Xe*  +  C0C1„  data.  The  fit  to  the  I_/I_  data 

u  D  C 

shown  in  Fig.  11,  though  not  as  good  as  for  the  other  two  reactions,  is  satis¬ 
factory.  The  distributions  based  on  simulation  of  the  Xe*  +  C0C17  spectra  are 


only  qualitatively  reliable  because  lack  of  oscillatory  structure  in  the 
spectra  prevents  reliable  assignment  of  the  vibrational  distributions.  Never¬ 


theless,  Model  C  does  give  acceptable  results  and  no  attempt  was'  made  to  further 

improve  the  model  for  the  low  v  regime. 

V.  RELAXATION  AND  TRANSFER  IN  OTHER  BATH  GASES 

V-A.  XeCl(high  v)  in  Helium  and  Neon 

In  this  section  we  extend  the  model  to  include  other  bath  gases.  For  the 

most  part  these  results  are  based  on  the  Xe*  +  reaction  and  will  be  strictly 

valid  only  for  XeCl  high  v.  Only  the  final  results  will  be  stated.  Examples 

of  the  results  of  the  simulation  and  modelling  results  for  these  gases,  such  as 

shown  in  Fig.  1-3  and  10-12  for  Ar,  can  be  found  elsewhere."^ 

Inspection  of  the  I  /I  ratios  (Fig.  4b)  for  XeCl  in  He  and  Ne  indicates 

d  L 

that  transfer  is  much  slower  in  these  bath  gases  than  in  Ar.  In  addition,  the 

XeCl(B,C)  vibrational  distributions  relax  much  slower  in  He  and  Ne  than 

in  Ar.  For  example  at  30  torr  of  He  the  distributiions  are  less  relaxed  than 

in  2  torr  of  Ar.  Before  discussing  the  He  and  Ne  results,  it  is  necessary  to 

£ 

make  a  few  comments  concerning  the  deduction  of  P„  .  For  the  Ar  case  sufficient 

changes  in  the  vibrational  distributions  occur  at  rather  low  pressures  so  that 

the  effects  of  just  a  few  collisions  with  the  bath  gas  can  be  observed  and  the 
£ 

elements  can  be  assigned  with  some  confidence,  i.e.  the  parameter  b  was 

determined  with  some  reliability.  For  He  or  Ne,  on  the  other  hand,  pressures  of 

10-20  torr  are  required  to  see  significant  changes  in  the  distributions  and 

small  changes  in  energy  with  large  cross-sections  are  difficult  to  distinguish 

from  larger  changes  with  smaller  cross-sections.  It  is  impossible  to 

£ 

simultaneously  determine  P  and  the  inelastic  collision  cross-section  with 

reliability.  For  C-B  transfer  this  situation  is  alleviated  because  the  I  /I 

B  C 

T 

ratio  identified  the  cross-section  (or  magnitude  of  k^)  and  the  model  for  P^ 
can  he  determined  more  reliably. 
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For  He  and  Ne  we  used  the  same  general  form  for  P. .  and  P. .  as  for  the  Ar 

ij  ij 

case  and  varied  the  parameters  until  a  successful  fit  to  the  data  was  found.  The 
best  trasnfer  model  for  He  was  with  C  =  0.2,  f  =  0.9  and  T.  =  0.15  with  T  -T,  the 

l  2^ 

O 

same  as  for  Model  C  and  a  collision  diameter  of  4  A.  This  gave  k  of  0.5-1. 8, 

1.8-2. 5,  and  2. 5-2. 5  x  10  ^cm^  molec  ^  sec  ^  for  the  low,  mid,  and  high  v  ranges, 
respectively.  For  vibrational  relaxation  a  range  of  parameters  could  be  used  with 
nearly  equal  success.  For  example,  with  R^-R^  the  same  as  f°r  m°del  C,  successful 
fits  were  obtained  with  b  =  0.1-0. 5  and  =  0.05-0.2.  However,  the  best  overall 
fit  was  judged  to  be  for  b  =  0.3  and  and  0.11  and  these  values  were  adopted 

for  the  best  fit  to  the  He  data.  These  gave  k  =  0.15-1.1,  1.1-2. 5,  and  2. 5-4.4  x 
10  ^  cm^  molec  ^  sec  ^  for  the  low,  mid,  and  high  v  ranges. 

The  neon  results  were  much  the  same  as  for  He.  Relaxation  again  cannot  be 
determined  accurately  and  several  models  could  fit  the  experimental  results 
satisfactorily.  The  best  model  (E)  used  Ri  =  0.3,  0.01,  1..0-0.3,  T^  =  0.3,  0.01,  0.7, 

-0.5  for  i  *  1  to  4  and  b  -  0.2,  c  -  0.15  and  d  =  0.87  with  a  collision  diameter 

0  -11  3 

of  5A.  These  parameters  gave  =  0.73-2.6,  2. 6-3. 5,  and  3- 5-3. 5  x  10  cm 

molec  ^  sec  ^  and  k  =  0.48-2.9,  2. 9-6. 2  and  6. 2-9. 5  10  cm^  molec  ^  sec  ^  for 

K 

the  low,  mid  and  high  v  ranges,  respectively. 

V-B.  XeCl(high  v)  in  Krypton,  Xenon  and  Nitrogen 

In  these  bath  gases  the  relaxation  and  transfer  rates  are  faster  than  in  Ar. 

The  simulated  distributions  could  be  represented  by  a  high  temperature  Boltzmann 
distribution  at  pressures  as  low  as  ^  5  torr.  The  problems  encountered  in  assigning 


the  relaxation  rate  constants  for  He  and  Ne  are  not  present  here  because  sig- 

T 

nificant  changes  in  the  distributions  are  observed  below  1  torr  and  both  P.. 

ij 

£ 

and  can  be  assigned  in  a  straight  forward  way.  For  Kr  we  used  the  same 
model  that  was  developed  for  the  argon  case.  The  best  parameters  are  b  =  0.5, 
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C  =  0.5,  f  =  0.8  and  =  0.2,  0.01,  1.0,  0.3  and  T\  =  1.1,  0.1,  0.7,  -0.5  for 

O 

i  *  1  co  4  and  with  a  collision  diameter  of  7A.  These  gave  k  =  0.6  -  2.7,  2.7  - 

5.5,  and  5.5  -  7.3  x  10  ^  cm"*  molec  ^  sec  ^  and  k  =  2.9  -11.0,  11.0  -  15.5, 

and  15.5  -  15.5  x  10  ^  cm^  molec  ^  sec  ^  for  the  low,  mid  and  high  v  ranges. 

The  Kr  transfer  rate  constants  are  only  slightly  larger  than  for  Ar  which  is 

consistent  with  the  I  /I  ratios  in  Fig.  4b.  However,  the  relaxation  rate  constants 

are  significantly  smaller  than  for  Ar  and,  in  fact,  are  even  less  than  for  N'e.  Yet 

the  observed  overall  rate  of  vibrational  relaxation  is  faster  in  Kr  than  in  Ar  because 

more  vibrational  energy  is  lost  per  collision.  Efforts  were  made  to  use  a  model 

with  smaller  loss  and  larger  k  to  fit  the  Kr  results  but  the  fit 

to  the  simulated  distributions  was  not  as  good  as  for  Model  F. 

Xenon  data  could  not  be  obtained  using  the  flowing  afterglow  or  resonance 

absorption  method  because  these  methods  fail  for  high  Xe  pressures.  Instead 

XeCl*  emission  was  obtained  in  the  presence  of  Xe  bath  gas  by  flowing  a  %  1% 

CI2  mixture  in  Xe  through  a  pulsed  discharge  and  observing  the  XeCl*  emission 

18 

directly  in  the  discharge;  the  details  will  be  reported  elsewhere.  This 
method  is  not  specific  in  the  excitation  mechanism;  i.e.  XeCl*  formation  may  result 
from  Xe(6s)  and  Xe  Rydberg  state  atoms  reacting  with  Cl.,  or  possibly  even 

4- 

ion-ion  combination.  Furthermore,  changing  the  Xe  pressure  may  alter  the  distribution 

13 

of  excitation  mechanisms.  Preliminary  interpretations  indicate  that  the  reaction 

of  CI2  with  Xe**  Rydberg  states  is  an  important  mechanism  for  the  discharge 

source  because  the  initial  vibrational  distribution  is  much  broader  than  the 
3 

one  from  the  Xe(  ,  reactions.  Although  a  detailed  analysis  of  the  vibrational 
relaxation  and  transfer  in  Xe  will  have  to  await  a  fuller  understanding  of  the 
excitation  mechanism,  we  estimate  from  the  I  /I  in  Fig.  4b  that  k  *■>  15-20  10  ^ 

D  C  I 

3-1-1 

cm  molec  sec  for  XeCl  (high  v)  +  Xe.  The  vibrational  relaxation 
resembles  the  Kr  case  with  (possibly)  even  more  vibrational  energy  lost  per 
collisions . 
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The  quality  of  the  experimental  data  for  XeCl*  in  was  poor  and  no 
detailed  calculations  were  performed.  However  the  appearance  of  the  spectra 


as  a  function  of  pressure  closely  resembles  the  Kr  bath  gas  case.  The 


I_/I_  ratio  (Fig.  4b)  is  also  similar  to  Kr.  Although  no  quantitative  results 

D  L 


can  be  given  it  appears  that  the  magnitude  of  the  rate  constants  and  the  form  of 


$  for  should  be  quite  similar  to  that  for  Kr. 


VI.  DISCUSSION 

VI. A.  Reliability  of  Rate  Constants  and  Models 

Beofre  interpreting  our  results  in  terms  of  molecular  interactions,  it  is 
necessary  to  assess  their  accuracy.  The  most  important  variable  affecting  the  magni¬ 
tude  of  the  rate  constants  is  the  lifetimes  used  in  the  calculations.  Since  these 
are  only  good  to  30%  and  since  the  population  distributions  are  uncertain  to 
+  10%,  the  magnitude  of  the  rate  constants  are  no  better  than  +  40%  on  an  absolute 


basis.  The  relative  values  should  be  better  than  this.  In  general  the  kT  should 


be  more  reliable  than  for  two  reasons:  1)  the  experimental  I^/I^  defines  the 


magnitude  of  k^  and  2)  the  simulated  distributions  of  the  B  state  are  more  reliable 


than  for  the  C  state.  Thus,  modeling  the  B  state  distribution  while  simultaneously 


fitting  the  IB/I_  ratio  in  the  low  pressure  regime  provides  a  reliable  estimate 
o  L 


for  kj,  (high  v) .  Of  course,  relaxation  and  transfer  are  not  independent  and, 


since  kj,  depends  on  vibrational  level,  the  rate  at  which  C  relaxes  will  have  an 


effect  on  the  rate  of  transfer.  This  was  noted  previously  with  Model  B,  which 
overestimated  k^  in  the  mid  v  range  causing  fast  vibrational  relxation  and  con¬ 


sequently  the  transfer  rate  decreased  too  quickly.  Nevertheless,  the  v  40% 


estimate  should  provide  a  reliable  limit  to  the  accuracy  of  our  results  for  k 


and  k^.  The  values  for  k^  (high  v)  deduced  from  modeling  agree  favorably  with 
the  estimates  from  a  two  state  model  (section  II)  affirming  that  the  high  v  level 
assignments  must  be  reasonably  reliable. 

With  the  above  comments  held  in  mind,  a  glance  at  Table  2  reveals  the  following 


-24- 


points  : 

(1)  The  k^  increase  in  the  order  He  <  Ne  <  Ar  <  Kr  <  ^  ''■>  Xe. 

(2)  The  k^  increase  in  the  same  order  except  for  Kr  (and  possibly  Xe). 

(3)  Both  k^  and  k^  depend  on  vibrational  level,  increasing  in  magnitude 
with  increasing  v. 

(4)  The  sum  of  k^  '+  k^,  is  nearly  gas  kinetic  in  magnitude  for  the  high  v. 

(5)  As  the  bath  gas  becomes  heavier,  k^  becomes  increasingly  important 

relative  to  k^  and  dominates  for  Kr  (ar.  ve?).  Fig.  13  shows  a  plot 

R  R 

of  z  P.,  (state  to  state  rate  constants)  vs.  AE  for  Model  A,  z  P..  and 
ij  iJ 

T 

z  P..  for  Models  B  and  C  for  v  =  108,  60  and  20.  For  v  =  108,  the 
13 

magnitude  of  z  P  is  similar  for  Models  A  and  B.  The  slight  difference 

in  shape  is  due  to  inclusion  of  the  1/AE  term  in  Model  A. 

Because  of  the  large  number  of  vibrational  levels  and  multiple  collisions  involved 

R  T 

in  a  steady-state  system,  unique  solutions  to  P..  and  P..  are  not  possible. 

T  R 

However,  the  models  we  developed  for  P. .  and  P..  are  useful  for  contrasting  the 

ij  ij 

relaxation  and  transfer  in  the  various  bath  gases.  Time  resolved  studies  would 

be  required  to  obtain  better  transition  probability  matrices. 

It  is  evident  from  Fig.  14  that  the  mean  energy  loss  per  transfer  collision 

T 

is  greater  than  for  relaxation;  however  the  variation  of  P,.  with  AE  is  nearlv 

13 

£ 

the  same  as  P  for  AE  >_  0.15  E  .  An  exponential  gap  model  could  not  be  used 
to  fit  transfer.  For  AE  <_  E^  it  was  necessary  to  reduce  the  transfer  in  this 
AE  range  so  that  the  B  vibrational  distribution  would  relax  fast  enough  and  still 


match  the  I  /I  ratio.  However,  the  bulk  of  the  transfer  takes  place  with  AE 
D  L 

similar  to  relaxation. 

R  T 

In  Fig.  15  the  difference  in  P..  and  P..  for  He  and  Kr  are  contrasted.  As 

13  ij 


can  be  seen,  the  results  are  quite  different  for  the  two  bath  gases. 


Both 


transfer  and  relaxation  in  Kr  involves  loss  of  large  amounts  of  vibrational  energy 

R  T 

and  significant  probability  for  large  quantum  transitions.  In  He  the  P  and  P_ 
are  grouped  more  sharply  around  the  initial  level.  Both  Xe  and  Ar  lie  between  the 


extremes  for  He  and  Kr. 

VI.  B.  Interpretations  of  Models. 

Before  discussion  of  rare  gas  -  XeCl  interaction  potentials  our  results 
will  be  compared  with  pertinent  data  in  the  literature.  Most  studies  of 
vibrational  relaxation  in  diatomic  molecules  have  focused  on  low  vibrational 
levels  where  the  spacing  between  levels,  IE,  is  greater  than  kT.  In  the  present 
case  the  inelastic  cross  sections  are  large  and  the  AE  <  kT.  For  this  reason 

4* 

we  compare  our  results  with  vibrational  deactivation  of  KBr  (+  represents 
vibrational  excitation)  formed  from  K  +  Br^-  The  relaxation  of  KBr  seems  an 
ideal  choice  for  several  reasons:  1)  reaction  of  Xe*  +  Cl^  is  in  many  ways^ 
similar  to  K  +  Br^  and  both  XeCl*and  KBr  have  <E^>  of  ^  40  kcal  mole  \  2) 
both  KBr  and  XeCl*  are  ionic  in  nature,  and  3)  an  extensive  amount  of  work  has 
been  reported  on  the  vibrational  deactivation  of  KBr.  Fisk  and  coworkers  studied 

KBr  relaxation  via  observation  of  inelastic  scattering  of  KBr  for  several  collision 

19  20 

partners.  In  addition  both  trajectory  studies  of  KBr  +  Ar  vibrational 

relaxation  and  an  information  theoretic  analysis  of  experimental  data  have  been 
21 

reported.  We  first  compare  our  results  with  their  data  by  the  information 

theoretic  approach.  Since  KBr  was  studied  by  beam  techniques,  the  experimental 

data  consists  of  the  flux  density  of  scattered  KBr  as  a  function  of  center  of 

19 

mass  scattering  angle  and  relative  velocity.  The  information  theoretic 

21 

analysis  of  Crim  and  Fisk  compared  the  measured  translational  energy 
distribution,  P(f^,),  where  f^  is  the  fraction  of  initial  energy  appearing  as 
product  translational  energy,  with  a  prior  of  reference  distribution,  P°(f^). 

The  appropriate  prior  for  diatomic-atomic  collisions  in  the  rigid  rotor-harmonic 
oscillator  limit  is 

P°(fT)  -  ATfT1/2  (l-fT)  (20) 

where  A_  is  a  normalization  constant. 
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Apparently  P(f^)  was  obtained  by  averaging  over  the  experimentally  observed  center 

21 

mass  angular  range  (v45  to  80°).  Since  there  is  a  significant  correlation 

between  the  scattering  angle  and  the  KBr  translational  velocity  distribution, 

these  "average"  experimental  distributions  may  not  be  good  representations  of 

the  total  distribution.  However  the  reported  distributions  are  clearly  two 

component  with  one  component  (the  larger)  centered  around  the  elastic  limit 

and  the  second  centered  around  a  large  inelastic  transfer.  These  two  components 

give  two  linear  regions  on  the  surprisal  plot.^2 

To  facilitate  the  comparison,  we  give  our  results  in  a  similar  manner. 

Because  the  initial  internal  energy  of  XeCl*  is  large  (E  v  40  kcal/mole 

mp 

for  Xe*  +  Cl^) ,  the  total  energy  of  the  rare  gas  +  XeCl*  collision  is 
approximately  fixed.  Furthermore  if  vibrational  to  translational  inelasticity 
is  the  dominant  deactivation  mechanism,  f  for  P..  is  given  by 

fT  =  E~j  +  ki- 
J 

The  distributions  P(fT)  and  P°(fT)  for  various  bath  gases  are  shown  in  Fig.  15 
for  XeCl  (v=100)  along  with  the  surprisal  1(f). 

I(fT)  =  -In  [P(fT)/P°(fT)]  (22) 

If  the  surprisal  is  linear,  the  departure  of  P(fT)  from  P°(f^,)  can  be 

described  by  a  parameter  1  where,  1  =  ttt^-T—  . 

d(fTJ 

If  '*0 ,  relaxation  takes  place  statistically  while  a  positive  or  negative  < 
indicates  relaxation  occurs  transfer  of  more  or  less  than  the  statistical  amount 
of  energy,  respectively.  As  can  be  seen  from  Fig.  15,  the  surprisals  are  linear 
for  He  and  Ne  but  deviate  from  linearity  for  Ar  and  Kr.  However,  a  linear  plot 
can  be  forced  and  the  1  values  given  in  Table  3  range  from  17.6  for  He  to  1.5 


for  Kr.  In  Fig.  13  it  was  shown  that  Model  A  and  B  predict  P..  which  deviate 

ij 

slightly  from  each  other.  This  deviation  is  largest  for  large  IE  transitions 


(i.e.  large  f^,)  and  consequently  the  surprisals  shown  in  Fig.  15  should  be  regarded 

as  uncertain  in  this  region.  Hence  the  linear  surprisal  parameters  were  obtained 

only  from  points  with  f ^  <0.7.  All  of  the  surprisals  show  a  marked  difference 

T 

(positive  deviation)  from  the  statistical  case.  A  similar  analysis  for  P _  gave 
the  same  qualitative  trends  for  the  bath  gases  but  with  A  reduced  by  a  factor  of 


v2  relative  to  P  .  The  large  positive  A  for  He  indicates  collisions  with  XeCl* 

are  nearly  elastic,  i.e.  little  energy  is  transferred  from  XeCl*  into  transition. 

Collisions  with  Kr  however  distribute  large  amounts  of  initial  vibrational  energy 

into  translation  and  bring  about  rapid  equilibration  of  translational  and 

internal  degrees  of  freedom.  Although  modeling  was  not  attempted  for  Xe  as  the 

bath  gas,  the  trends  evident  in  Table  3  would  predict  that  Xe  -  XeCl*  collisions 

bring  about  a  statistical  or  near  statistical  exchange  of  energy. 

The  results  of  the  surprisal  analysis  are  contrasted  in  Table  3  with  the 

21 

results  of  vibrational  relaxation  of  KBr.  As  already  noted  Crim  and  Fisk 

found  two  linear  regions  in  their  surprisa]*  plot  for  inelastic  scattering  of  KBr. 

The  component  at  low  fT  gave  a  large  A  consistent  with  nearly  elastic  scattering 

and  the  f^,  >  0.2  component  gave  A  "u  0  which  corresponds  to  nearly  statistical 

redistribution  of  energy.  The  I(f^)  in  Figure  15  for  Ar  and  Kr  also  could  be 

approximated  by  a  two  component  distribution  similar  to  that  found  for  KBr. 

The  distributions  of  Crim  and  Fisk  strongly  resemble  the  sum  of  an  exponential 
-AE  / 

P .  .  a  e  ij  b  and  linear,  P..  *  aAE..,  distribution.  However,  theirs  is  much 
ij  ij  U 

more  extreme  form  than  our  model  A.  We  must  conclude  that,  although  the  two 
transition  probability  models  are  qualitatively  similar,  they  differ  significantly 
at  the  quantitative  level.  This  disagreement  could  be  a  consequence  of  the  lack 
of  inclusion  of  all  scattering  angles  in  the  beam  studies,  the  difficulty  in 
obtaining  single  collision  transition  probabilities  from  our  steady-state 
experiments,  or  a  real  difference  as  a  consequence  of  the  different  intermoiecular 
potentials  (see  next  section)  for  KBr  and  XeCl  +  rare  gas. 


Included  in  Table  3  are  Che  mean  energy  loss  collisions  for  XeCl  relaxation 
and  transfer  and  for  KBr  relaxation.  For  Ar  and  Kr ,  the  <AE>  for  XeCl  is  the 
order  observed  for  KBr  deactivation.  However,  Me  appears  to  deactivate  KBr 
much  more  efficiently  than  XeCl*.  In  fact  Me  appears  to  deactivate  KBr 
more  efficiently  than  Ar,  which  is  contrary  to  what  one  would  expect.  Simple  theory 
as  well  as  experimental  evidence  indicate  that  for  efficient  coupling,  the 
collision  time  and  the  vibrational  period  of  the  oscillator  should  be  closely 
matched.  For  light  molecules  in  low  vibrational  levels  with  large  vibrational 
frequencies,  this  results  in  lighter  gases,  eg.  He,  being  more  efficient  for 
vibrational  relaxation.  For  KBr  or  XeCl*  in  high  vibrational  levels,  the 
oscillator  frequency  is  small  and  slower  collisions  may  be  more  efficient  for 

deactivation.  The  mean  collision  time  x^  was  estimated  for  an  interaction  length 

0  1/2 

of  'u  2  A  and  the  mean  velocity  at  300K,  v  =  (8kT/FIp)  .  For  XeCl* (v=100) the 
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vibrational  period  xv  is  ^  4.5  x  10  sec,  and  xc/x^,  is  'u  0.3,  0.5,  0.7,  1.0, 

and  1.1  for  He,.Ne,  Ar,  Kr,  and  Xe, respectively.  The  match  of  x^  and  xv  is  best 

for  Kr  and  Xe  which  is  consistent  with  larger  <AE>  results  in  Table  3.  Of 

course,  the  interaction  potentials  also  are  more  attractive  for  the  heavier  rare  gases. 

In  an  effort  to  determine  the  features  of  the  interaction  potential  that  wer« 

important  to  KBr  vibrational  deactivation,  Matzen  and  Fisk  performed  a  classical 

20 

trajectory  study  for  the  KBr  +  Ar  system  using  three  different  potentials, 

Surface  I,  the  best  surface  was  constructed  empirically  from  the  two-body 
interactions  of  K  ,  Br  and  Ar  and  an  Ar-KBr  interaction  term.  Surface  II  was 
like  Surface  I  but  ignored  all  attractive  interactions.  Surface  III  was  similar 
to  II  but  substituted  a  spherically  symmetric  attractive  part  for  the  Ar  +  KBr 
interaction.  The  major  differences  in  the  potentials  were  the  attractive  inter¬ 
actions.  Surface  II  was  purely  repulsive  while  the  attractive  features  of  I 
and  II  were  of  the  same  strength  but  of  different  shape.  On  I  the  incoming  Ar 
is  attracted  preferentially  to  the  K+  end  of  KBr  while  in  III  the  interaction  is 


symmetric . 


The  results  and  conclusions  of  the  trajectory  study  which  are  pertinent  to 
XeCl*  relaxation  can  be  summarized  as  follows: 

i)  Both  the  shape  and  strength  of  the  attractive  interactions  were 
important.  Surface  I  gave  the  best  agreement  with  experimental 
results  and  the  calculated  distribution  did  have  a  two  component 
nature. 

ii)  The  strength  of  the  attractive  interaction  determined  the  magnitude 
of  the  inelastic  cross-section  and  both  I  and  III  gave  higher 
calculated  cross  sections  than  II. 

iii)  The  shape  of  the  attractive  potential  was  important  in  determining 

the  average  amount  of  energy  transferred.  For  KBr  +  Ar,  collisions 
with  the  K+  end  resulted  in  larger  inelasticities  than  collisions  with 
the  Br  end  and  surface  I, which  attracted  the  incoming  Ar  to  the  K+  end 
of  KBr?gave  larger  inelasticities. 

iv)  The  high  energy  component  was  associated  with  strong  repulsive  interactions 
Although  attractive  interactions  were  important  in  bringing  the  collision 
partners  together,  large  inelasticity  was  shown  only  for  trajectories 
which  sampled  the  strongly  repulsive  regions  of  the  potential  energy 
surface. 

The  XeCl  situation  will  be  more  complicated  than  for  X5r  -  Ar  because  of  multiple 
surfaces  with  differing  degrees  of  attractive  and  repulsive  character  and  because 
of  the  coupled  vibrational  relaxation  and  electronic  state  transfer  channels. 
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Although  accurate  potentials  are  not  known  for  XeCl  -  Rg ,  qualitative 

23  ?4 

features  can  be  estimated  from  ab-initio  and  diatomic-in-molecules 
calculations  which  have  been  reported  for  some  homonuclear  trimer  systems, 

Rg£X  (Rg  =  Ar,  Kr  and  X  =  F,  Cl).  The  three  lowest  states  of  Rg2X  are  covalent 
and  repulsive  and  correlate  to  RgX(A  or  X)  and  Rg;  these  states  are  of  no 

importance  here.  The  next  three  higher  states  are  ionic  in  nature  and  bound 

+  -  ? 
relative  to  Rg^  +  X  .  Only  the  lowest  ionic  state  (2  in  C^v  geometry) 

"k 

though,  is  bound  with  respect  to  RgX  +  Rg.  The  most  stable  geometry  of  this  lowest 

ionic  state  is  an  isoceles  triangle.  For  the  mixed  rare  gas  trimers,  such  as 

RgXeCl,  the  above  potentials  will  be  altered  somewhat.  First,  the  isoceles 

triangle  geometry  will  be  distorted  with  the  Rg  atom  drawn  closer  to  the  Xe  end 

of  XeCl.  Secondly,  RgXeCl  trimers  would  be  expected  to  be  much  less  strongly 

+  + 

bound  than  Xe2Cl,  paralleling  the  decreasing  binding  energy  of  RgXe  vs  Xe2  • 

*  *  25  26 

In  fact,  emission  from  both  Xe2Cl  and  KrXeCl  have  been  observed.  ’  Of  the 
*  * 

lighter  RgXeCl  ,  only  ArXeCl  would  be  (possibly)  weakly  bound  by  more  than 

van  der  Waals  forces.  According  to  the  calculations,  the  bound  Rg2+X  state 

2  A 

correlates  with  RgX(C)  and  the  next  higher  ionic  (3  B„  in  C~  or  A'  in  C 

l  Zv  s 

geometry)  state  with  RgX(B).  The  approach  of  Rg  to  XeCl(C)  in  Cg  geometry 

gives  an  A'  and  an  A"  state  (neglecting  spin-orbit  interaction)  depending  on 

whether  the  sideways  approach  of  Rg  to  XeCl  is  head  on  or  perpendicular  to  the 

singly  occupied  II  orbital.  The  approach  of  Rg  to  XeCl(B)  gives  only  an  A' 

state  as  shown  in  Fig.  16.  If  the  C  state  is  higher  in  energy  than  the  B  state, 

this  results  in  a  curve  crossing  which  can  lead  to  collisional  coupling  of  RgX(B,C). 

Since  XeCl(C)  is  actually  lower  in  energy  than  XeCl(B),  the  Rg  +  XeCl(C)  potential  shou 

*  2 

correlate  adiabaticallv  to  the  lowest  ionic  state  of  RgXeCl  (3“A'  in  C  geometry) 

s 

as  shown  in  Fig.  16a.  The  Rg  +  XeCl(B)  potential  is  repulsive  and  the  interaction 
of  the  two  A'  states  should  be  small  and  collisional  mixing  of  XeCl*(B,C)  should 
be  slow, which  is  contrary  to  observation.  However, the  ab  initio  calculations  show 


that  the  collinear  approach  of  Rg  to  RgX*  is  also  attractive,  although  the 

binding  is  not  as  large  as  for  the  approach,  as  shown  in  Fig.  16.  In  this 

case  approach  of  Rg  to  XeCl(B)  is  attractive,  while  Rg  +  XeCl(C)  is  repulsive. 

However  the  states  are  of  different  symmetry  and  coupling  of  XeCl(B,C)  via  this 

approach  would  also  be  weak  unless  spin-orbit  interaction  is  included.  For  such 

heavy  atoms  as  Kr  or  Xe  this  could  be  of  importance  and  lead  to  mixing  of  XeCl(B,C) 

Even  though  spin  orbits  effects  allow  mixing  of  XeCl(B,C)  the  efficiency  of  this 

transfer  is  still  surprising.  Furthermore,  the  rapid  B-C  transfer  holds  for  all 
3  4 

RgX  molecules  ’  regardless  of  the  B  and  C  state  ordering  and  the  magnitude  of 
the  spin-orbit  interactions .  Thus,  the  asymptotic  properties  of  Fig.  16  are  not 
critical  for  highly  vibrationally  excited  molecules  (they  may  well  be  of  importance 
for  v=0-5  levels  when  the  transfer  rates  are  1-2  orders  of  magnitude  smaller). 

Fig.  16  also  portrays  the  behavior  of  the  adiabatic  potentials  as  a  function  of 
ot,  the  angle  between  the  Rg  -  Rg  -  X  bonds.  As  a  increases  from  v  60°  to  the 
linear  configuration,  180°,  the  3  A'  and  4  A'  give  rise  to  an  avoided  crossing. 

This  crossing  will  exist  for  all  RgX*  +  Rg  and  presumably  is  responsible  for 
electronic  state  transfer.  Apparently  the  large  amplitude  motion  on  the  potential 
curves  by  highly  vibrationally  excited  RgX*  will  lead  to  frequent  sampling  of 
this  region  where  the  potentials  interact.  Furthermore,  collisions  of  both 
RgX(B)  and  RgX(C)  can  sample  the  lowest  potential  depending  on  the  approach  angle 
of  Rg. 

The  KBr  -  Ar  trajectory  study  found  that  the  strength  of  the  attractive 
interaction  determined  the  magnitude  of  the  inelastic  cross  sections.  This 
behavior  is  generally  followed  here  for  He  -  Ar .  However,  for  Kr  the  relaxation 
rate  constants  were  much  slower  with  C-B  transfer  lominating.  Only  for  Kr  (and 
Xe)  would  the  3  A'  potenital  be  significantly  bound  and  transfer  may  dominate 
over  relaxation  for  this  reason. 


Simple  theory  would  predict  that  Ar  striking  the  K  end  of  KBr  would 

result  in  larger  inelasticities  than  Ar  striking  the  Br  end  because  K+  is  much 

—  4- 

lighter  thanBr  .  In  XeCl*  the  attraction  is  to  the  Xe  end  rather  than  the 

lighter  Cl  and  our  results  show  somewhat  smaller  inelasticity  for  Ar  or  Ne  than 

in  KBr.  Because  of  the  stronger  attractive  interactions  and  presumably  softer 

repulsions  for  XeCl*  the  present  results  suggest  that  the  attractive  interactions 

27 

play  a  more  dominant  role  in  XeCl*  +  Rg  that  KBr  +  Rg.  Greene  et  al,  in 

studying  the  inelastic  scattering  of  Csl  by  Ar  and  Xe  suggest  that  short  lived 

complex  may  be  formed  at  low  collisions  energies.  Their  results  fit  statistical 

27 

partitioning  of  the  total  energy.  Greene  et  al  also  suggest  that  in  the  inter¬ 
action  of  a  highly  excited  ionic  molecule,  such  as  KBr  or  Csl  with  Ar ,  that 
during  the  bond  extension  primve  of  the  vibration  the  interaction  may  resemble 
Ar  interacting  with  K+.  Such  an  interaction  is  much  stronger  than  with  the 
molecule  leading  to  stronger  attractive  interaction  during  the  collision.  For 
XeCl*  +  Rg  this  effect  will  be  enhanced  because  RgXe+  is  chemically  bound  (at 
least  for  Rg  =  Kr,Xe)  whereas  in  ArK+  the  attraction  is  only  electrostatic  in 
nature.  We  see  some  evidence  for  this  because  of  the  dependence  of  the  rate 
constants  on  vibrational  excitation.  For  the  high  v  the  attractive  interaction 
is  stronger  and  leads  to  larger  inelastic  cross  sections  as  opposed  to  the  low 
v  levels  where  the  interaction  is  more  like  the  weaker  Rg-XeCl  interaction. 

VI.  Conclusions 


From  computer  simulation  of  the  XeCl(B-X)  and  XeCl(C-A)  bound-free  emission 
spectra  in  He,  Ne,  Ar,  and  Kr  bath  gases,  the  vibrational  and  electronic  state 
distributions  were  assigned  for  a  range  of  pressures.  S tate-to-state  rate 
constants  for  vibrational  relaxation  and  electronic  state  transfer 
have  been  assigned  for  XeCl*(B,C)  +  Rg  collisions  by  master  equation 
calculations  fitting  the  variation  of  the  vibrational  distributions  with  Rg 


pressure.  The  electronic  state  transfer  rate  constants  for  the 
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cm  molec  sec  for  the  Kr-He  series;  the 


range  varied  from  15-2.5  x  10 


the  vibrational  relaxation  rate  constants  are  comparable  with  the  exception  of 
Kr  which  is  smaller  than  for  Ar.  The  magnitude  of  both  the  relaxation  and  transfer 
rate  constants  decline  by  an  order  of  magnitude  from  v=100  to  v  =  10.  The  model 
for  state-to-state  vibrational  relaxation  is  adequately  represented  by  the 
exponential  gap  law.  For  electonric  state  transfer  the  change  of  vibrational 
energy  is  larger  than  for  the  collisions  giving  vibrational  relaxation  without 
change  in  electronic  state,  and  the  model  is  Poisson  in  nature.  The  larger  loss 
of  vibrational  energy  associated  with  electronic  state  changing  collisions  is 
explained  by  the  attractive  potential  involved  in  interconversion  of  the  XeCl(B) 
and  XeCl(C)  states. 
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Table  1.  XeCl  Energy  Levels  and  Radiative  Lifetimes 


Vibrational 

C  state 

B  state 

Level 

Energy  (cm 

i  (nsec) 

Energy  (cm 

'  (nsec 

0 

139.42 

120.0 

369.42 

11.0 

4 

876.08 

127.6 

1136.05 

11.08 

8 

1590.86 

136.4 

1882.33 

11.88 

12 

2284.25 

137.2 

2608.63 

12.29 

16 

2956.77 

142.8 

3315.38 

12.64 

20 

3608.94 

146.9 

4002.98 

12.53 

24 

4241.29 

152.3 

4671.84 

12.35 

28 

4854.33 

174.1 

5322.39 

13.43 

32 

5448.6 

182.1 

5955.05 

14.10 

36 

6024.61 

195.3 

6570.25 

14.50 

40 

6582.89 

195.5 

7168.42 

14.84 

44 

7123.96 

210.3 

7750.00 

16.12 

48 

7648.33 

224.6 

8315.41 

16.38 

52 

8156.52 

230.6 

8865.10 

17.25 

56 

8649.03 

245.0 

9399.49 

18.69 

60 

9126.35 

256.4 

9919.03 

19.33 

64 

9588.98 

265.0 

10424.17 

20.15 

68 

10037.4 

275.2 

10915.27 

21.35 

72 

10472.1 

279.1 

11392.77 

22.42 

76 

10883.4 

270.2 

11897.07 

23.88 

80 

11302.0 

296.2 

12308.67 

24.78 

84 

11698.1 

298.2 

12747.97 

26.04 

88 

12082.3 

308.3 

13175.27 

27.52 

Table  1.,  Continued 


t-  iwiiax 

Level 

c  state 

B  state 

Energy  (cm  ) 

i  (nsec) 

Energy  (cm  *) 

(nsec 

92 

12454.9 

318.1 

13590.97 

28.98 

96 

12815.3 

325.6 

13995.47 

30.21 

100 

13167 

337.7 

14389.17 

31.77 

104 

13507.3 

343.3 

14772.37 

33.21 

108 

13837.6 

349.1 

15145.147 

35.14 

112 

14158.1 

352.8 

15508.67 

37.16 

116 

14469.3 

357.9 

15862.27 

39.03 

120 

14771.5 

375.1 

16206.67 

40.91 

124 

15065 

398.5 

16541.97 

128 

15627.1 

433.7 

17186.47 

136 

15896.2 

438.5 

17496.07 

Rate  Constant  Summary 
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FIGURE  CAPTIONS 


Figure  1.  Comparison  of  experimental  (  -  )  and  simulated  ( - )  spectra 

of  the  B-X  and  C-A  emission  of  XeCl  resulting  from  the  reaction 
of  Xe*  +  Cl0  in  Ar.  In  each  case  the  B-X  spectra  have  been 
attenuated  so  that  the  most  intense  peak  of  the  C-A  and  B-X 
emission  have  the  same  height. 

Figure  2.  Same  as  Figure  1  but  from  Xe*  +  CCl^.  In  the  240-260  nm  region, 

the  simulated  B-X  spectra  are  above  the  experimental  spectra 
because  a  flat  potential  was  used  for  the  XeCl(X).  The  use 
of  a  more  realistic  potential  would  reduce  the  intensity  by 
shifting  the  transition  occuring  at  the  inner  turning  point  of 
the  upper  state  to  longer  wavelengths. 

Figure  3.  Simulated  ( - )  and  experimental  ( - )  spectra  from  Xe*  +  COC^. 

The  lack  of  oscillatory  structure,  even  at  low  pressure  prevents 
a  reliable  viorational  distribution  from  being  assigned. 

Figure  4.  Comparison  of  the  I^/I^  intensity  ratios  as  a  function  of  pressure 

for  a)  The  CCl^  and  COC^  donors  in  argon  and  b)  for  Xe*  +  Cl 

in  the  various  bath  gases.  The  dashed  line  in  b  serves  only  as  an 
aid  in  identifying  the  argon  points. 

Figure  5.  Steady-state  XeCl(B,C)  vibrational  distribution  obtained  from  simu 

lation  of  Xe*  +  spectra  at  various  Ar  pressures.  Because  of 

different  w  and  u>  X  values,  the  v  levels  for  the  two  states 
e  e  e  max 

are  different;  even  if  the  vibrational  energy  is  the  same,  see 
Table  1. 

Figure  6.  Steady-state  XeCl(B,C)  vibrational  distribution  obtained  from 

simulation  of  the  Xe*  +  CCl^  spectra  at  various  pressures  of  Ar. 

Figure  7.  Steady-state  XeCl(B,C)  vibrational  distributions  from  Xe*  +  C0C1., 

at  0.1-0. 2  and  3-6  torr.  Only  the  shapes  of  the  distributions 
have  significance;  note  the  similarity  to  3oltzmann  distributions. 
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Figure  11. 
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Figure  14. 


a)  Comparison  of  XeCl(B-X)  experimental  spectrum  with  two 

simulated  spectra  ( -  and  ••••)»  which  give  satisfactory  fits. 

The  simulated  spectra  for  the  Cl^  case  are  nearly  superimposed. 

b)  Corresponding  vibrational  distributions. 

Comparison  of  simulated  distributions  ( - )  and  calculated  distri¬ 
butions  from  Model  A  ( - )  for  the  vibrational  relaxation  of 

XeCl(C)  in  Ar. 

Comparison  of  simulated  ( - )  and  calculated  distributions  from 

Model  B  ( - )  of  XeCl(B)  and  XeCl(C)  from  Xe*  +  Cl  in  Ar.  Model 

B  includes  both  relaxation  and  transfer  processes.  The  apparent 
populations  in  the  lowest  two  levels,  which  represent  the  real 
v=0-8  levels  for  the  higher  pressures  can  not  be  interpreted 
literally,  see  text. 

Comparison  of  experimental  and  calculated  I„/I„  ratios  in  Ar.  For 
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clarity  the  I^/I^  ratios  have  been  moved  up  4  and  2  units  for  Cl„ 
and  CCl^,  respectively.  The  calculated  results  are  the  dashed  lines 
(best  fit)  or  dotted  line.  The  best  fit  results  are  from  Model  B 
(C^)  and  Model  C  for  CCl^  and  COC^*  The  dotted  line  is  a  Model  B 
calculation  for  CCl^;  the  poor  fit  necessitated  an  improvement  and 
led  to  Model  C. 

Comparison  of  Model  B  (•••)  and  Model  C  ( - )  results  for  Xe*  +  CC1 

data.  The  poor  fit  for  Model  B  required  improvement  and  Model  C 

was  developed  to  fit  the  Xe*  +  CCl^  data  (see  text). 

Comparison  of  relative  z-P„  f°r  Model  A  and  B  for  starting  v=20 

40,  and  108.  •  =  zP^  (Model  A),  ■  -  Z  (Model  B)  ,  and 

★  =  zP^  (Model  B)  . 
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Comparison  of  P^  (upper  panel)  and  P  (lower)  for  He  (Model  D) 

•and  Kr  (Model  F)  as  bath  gases.  To  emphasize  the  differences 
between  the  two  models,  the  maximum  P  have  been  noraaiized  to  10 


Figure  15. 


Comparison  of  translational  energy  distributions  P(f  )  ( - ) 

(lower)  and  the  corresponding  surprisals  (upper)  ,  calculated 
from  Eqn.  22  for  the  various  bath  gases.  Because  the  calculations 
included  only  every  fourth  level  the  origin  of  f  is  0.05. 

Figure  16.  Schematic  RgRgX*  potentials  as  a  function  of  a,  the  Rg-Rg-X 

band  angle  and  the  separated  Rg  +  RgX*  states  to  which  they 
correlate  .  The  potential  is  drawn  for  a  bound  trimer,  i.e. 
KrXeCl*  or  Xe2Cl*.  The  C9^  geometry  is  realized  only  for 
Xe2Cl*.  For  the  lighter  rare  gases  all  the  curves  become 
less  attractive  and  finally  strictly  repulsive  in  the  series 


Ar  <  Ne  <  He. 
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